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SUMMARY 


The  transverse  pumping  of  a  cylindrical  laser  rod  such  as  YAG 
by  a  set  of  diode  laser  arrays  equally  spaced  in  angle  around  the  rod 
has  been  analyzed  using  both  geometrical  and  wave  optics.  Between  the 
pump  sources  and  the  laser  rod  there  may  be  an  arbitrary  number  of 
concentric  media,  transparent  to  the  pump  light;  and  in  addition 
modification  of  the  diode  laser  radiation  patterns  will  require  other 
optical  components.  The  analysis  accommodates  both  of  these 
complications,  even  when  the  optics  in  front  of  the  diodes  do  not  image 
point  sources  to  points.  It  is  sufficient  to  be  able  to  describe  the 
intensity  and  ray  directions  in  finite  apertures  concentric  with  the 
rod. 


One  result  of  the  cylindrical  geometry  is  the  formation  of 
caustics  which  have  infinite  intensity  in  the  geometrical  optics 
approximation;  and  it  is  necessary  to  develop  a  wave  optics  theory  to 
predict  the  intensity  in  the  neighborhood  of  the  caustic  edge  of  the 
illuminated  region  inside  the  laser  rod.  The  wave  theory  first  uses  the 
geometrical  optics  method  to  calculate  the  amplitude  and  phase  just 
inside  the  rod's  surface,  which  serves  then  as  an  aperture  distribution 
function  describing  the  source  of  Huygens  wavelets  in  the  computation  of 
a  diffraction  integral An  approximation  that  assumes  constant  ampli¬ 
tude  and  cubic  phase  distribution  along  the  source  surface  yields  an 
Airy  function  description  if  the  caustic  intensities  that  is  useful  for 
relating  the  problem  to  other  known  physical  phenomena,  but  over¬ 
estimates  the  intensity  peak  significantly  in  this  case  when  compared 
with  the  results  of  numerical  computation  of  the  integral.  Computation 
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of  the  diffraction  integral  is  a  more  convenient  method  of  estimating 
the  peak  intensity  near  the  caustic. 

A  problem  with  transverse  pumping  is  that  if  the  absorption  is 
small  enough  to  pump  the  axial  region  of  the  rod  effectively,  a 
significant  fraction  of  the  pump  light  is  transmitted  through  the  rod. 
This  problem  can  be  overcome  by  reflecting  the  transmitted  light  back  to 
the  rod.  We  assume  that  the  reflecting  surface  is  at  the  same  radius  as 
the  source  images  or  source  apertures.  The  analysis  shows  that  this 
measure  leads  to  a  high  pumping  efficiency,  defined  by  the  fraction  of 
pump  light  absorbed  in  the  laser  rod.  All  of  our  computations  of 
pumping  intensity  have  assumed  a  Gaussian  profile  for  the  source 
radiation  pattern.  All  pump  light  outside  the  limiting  rays,  which 
graze  the  laser  rod,  is  wasted.  This  is  the  principal  efficiency 
limiting  mechanism  according  to  the  present  analysis,  and  this  points  to 
the  need  for  mirrors  or  specially  formed  lenses  near  the  sources  to 
redirect  the  radiation. 

The  analysis  has  been  simplified  in  many  portions,  and  is 
presented  in  a  form  that  makes  it  easy  to  program  on  a  computer .  The 
main  parts  that  we  did  not  program  are  those  that  would  account  for 
longitudinal  components  in  the  pump  intensity ,  although  we  have 
indicated  how  this  could  be  done.  The  two  dimensional  theory  will  be 
adequate  for  most  puposes,  specifically  in  computations  that  aim  at 
optimizing  the  design  of  a  laser  by  maximizing  some  linear  combination 
of  pumping  efficiency  and  deviation  from  a  desired  pumping  profile. 
Programs  have  been  written  for  the  IBM  PC  for  computing  intensity 
distributions  using  both  geometrical  and  wave  optics,  as  well  as  for 
investigating  numerically  some  of  the  relations  amongst  design  para¬ 
meters. 
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SECTION  1 


TRANSVERSE  PUMPING  -  GEOMETRICAL  THEORY 


This  section  treats  the  transverse  pumping  of  a  cylindrical 
laser  rod  by  diode  lasers.  The  aim  is  to  develop  appropriate  mathemat¬ 
ical  formulae  to  calculate  the  pump  energy  distribution  in  the  rod, 
determine  what  factors  influence  the  uniformity  and  efficiency  of  pump 
energy  deposition,  and  provide  some  logical  basis  for  the  choice  of 
dimensions,  number  of  arrays  on  the  circle,  etc.  The  simplest  approxi¬ 
mation  is  to  treat  the  problem  in  two  dimensions  (transverse  to  the  rod 
axis)  and  this  is  done  first;  but  the  problem  is  really  three-dimen¬ 
sional,  and  the  relevant  equations  are  developed  for  this  case  also. 

The  geometry. 

Figure  1  shows  a  section  of  the  arrangement  of  diode  arrays 
pumping  a  rod.  The  arrays  are  equally  spaced  on  the  circle.  Each  array 
has  an  associated  antenna  pattern.  A  phased  array  would  have  a  two-lobe 
pattern,  and  the  arrays  would  probably  be  arranged  so  that  the  lobe 
separation  is  in  a  direction  perpendicular  to  the  figure.  Each  array 
may  have  a  lens  in  front  of  it.  The  lenses  have  three  effects;  addi¬ 
tional  attentuation  due  to  surface  reflections ,  transformation  of  the 
antenna  pattern,  and  replacement  of  array  positions  by  image  posi¬ 
tions.  The  lens  arrangement  may  produce  real  or  virtual  images.  Figure 
1  is  still  good.  The  antenna  patterns  and  array  positions  are  simply 
those  of  the  images.  All  formulae  developed  initially  are  for  the  two- 
dimensional  approximation,  i.e.  all  rays  are  assumed  to  lie  in  the  plane 
of  the  figure. 
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Snell's  law. 


In  two-dimensional  coordinates  (r,B)  ,  when  the  index  n  is 
dependent  only  on  r,  Snell  's  law  is 

n  r  sinB  =  const. 

where  8  is  the  angle  between  the  ray  and  the  radius.  A  point  inside  the 

laser  rod  is  denoted  by  (a,  a),  and  is  referred  to  as  the  field  point. 

Figure  2  shows  how  the  field  point  is  illuminated  by  a  ray  at  angle  4> 

from  the  array  located  at  (R,Q)  .  Concentric  media  with  radii  ■ 

.  .  .rM  and  refractive  indices  n n^,  .  .  .nH  are  shown,  with  the  ray 

having  angles  of  incidence  and  refraction  8'.  and  8^.  at  the  point 

fr.,9.).  The  coordinates  (r  ,  9  )  and  (r„,  9.J  are  to  be  identified  with 

j  j  o  o  MM 

(a,  a)  and  (R,Q  )  respectively  when  the  index  j  takes  on  its  extreme 

M 

values  0  and  M.  The  angle  0  is  the  position  of  the  diode  laser  array, 

and  will  be  denoted  by  0 ^  where  i  =  1,  2,  .  .  .  .N  for  N  arrays.  The 

angle  4>  at  which  a  ray  leaves  the  diode  array  is  identified 

with  8  ;  and  finally,  we  also  use  n  interchangeably  with  n„.  The 
M  n 

closest  approach  of  the  ray  to  the  center  of  the  rod  is  denoted  by  the 

distance  p.  Snell's  law  applied  to  Figure  2  gives: 

n^p  =  n  R  sin4> 

=  n  .r  .  sinQ  . 

J  J  J 


-  n  .  ,  r  .  sin B'.  , ,  , 

j+1  j  j  (1) 

These  relations  give  the  various  angles  in  terms  of  $  .  In  general,  a 
ray  will  cut  the  circle  of  radius  a  at  two  points:  (a,  a^)  and  ( a ,0^), 
as  shown  above  in  Figure  3.  Identifying  the  ray  <t>  that  goes  through  the 
field  point  (a, a)  involves  expressing 
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Figure  3.  Illustration  of  how  the  ray  Intersects  the  circle  of  radius  a 
at  two  points  with  polar  angles  of  and  a^. 
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and  in  te  ms  of  <J>  .  The  following  relations  between  the  8'.,  8'.  and 

9  .  : 

J 


6'.  ,-8.  =  9.  ,  -  9. 
J-i  J  J-l  J 


(2) 


which  can  be  summed  over  j  to  give 


M 

9  -9=1  (B'  .  -  8  J 

o  M  .=1  j-l  j 


(3) 


where  6'  is  seen  to  be  sin  ^  (p/a)  ,  which  by  Equation  1  is  the  same  as 
sin  *  (^—  sing) .  Note  also  from  Figure  2  that 

n  ia  3 


a2  =  *2  cos  1(p/a) 


(4) 


Combining  Equations  1,  3  and  4  gives: 


.  *1  .  .  -1 ,  nR 

9  +  I  { sm  (—— - 

n  .  n  .r  .  , 

j=i  j  j-i 


sing)  -  sin  * (  nR  sing)} 
n  .r  . 

J  J 


(5a) 


and 


_  .  -1 ,nR  ,  . 

a  „  =  a ,  +  n  -  2  sm  ( -  cos g) 

2  1  n  a 


(5b) 


Notice  that  the  last  term  in  Eq.  5a  reduces  (at  j  -  H)  to  <t>  . 

As  $  increases,  8'  eventually  becomes  i \/2  which  determines  the 
limiting  ray  beyond  which  no  light  enters  the  rod.  From  Equation  1, 
this  condition  is  expressed  as 


TTC-1 529-R 


6 


(6a) 


assuming  that  n j  is  greater  than  n2.  Light  emitted  at  angles  exceeding 

<t>  at  the  diode  array  is  wasted  and  this  is  one  reason  for  considering 
m 

the  use  of  lenses.  (In  the  present  analysis,  a  Gaussian  distribution 
2 

exp  -  has  been  assumed.)  .A  lens  at  a  distance  u  from  the  point 

source  (section  through  diode  array)  produces  an  image  at  v  from  the 
lens  so  that  the  antenna  pattern  is  narrowed  by  the  factor  u/v  while  R 
increases  by  -(u  +  v) .  (v  is  positive  for  a  real  image.)  The  focal 
length  and  position  of  the  lens  would  be  specified  so  as  to  place  the 
image  at  the  desired  R  and  make  4>e  equal  to  the  desired  value. 

Equations  1-5  above  are  the  basic  ray  trac'ng  relations  in 
two  dimensions .  Eq.  5  can  be  used  to  find  the  polar  angle  a  at  which  a 
given  ray  from  the  source  intersects  the  circle  r  =  a  inside  the  rod. 
The  number  of  intersections  is  2,  1  or  0  depending  on  whether  a  is 

greater  or  smaller  than  pm. 

Figure  4  is  a  plot  of  Eqs .  5  for  a  specific  concentric 

geometry  chosen  for  purposes  of  illustration.  Here  <(>  is  6.07  degrees, 

nt 

for  the  critical  ray.  Curves  are  drawn  for  six  separate  values  of 
radius  a,  with  and  being  the  lower  and  upper  branches  respec¬ 
tively.  The  branches  are  separate  when  a  is  greater  than  (n2/n q)r^  and 
form  a  single  curve  for  smaller  values.  Curve  4  has  a  equal  to 
( n 2/n ^) r 2  which  corresponds  to  the  case  in  which  the  critical  ray  is 
tangential  to  the  circle.  The  remaining  curves  from  left  to  right  have 
a/r2  equal  to  0.2,  0.4,  0.6,  0.8,  1.0  respectively .  The  curvature  of 
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the  upper  branch  ( )  changes  sign  as  a  function  of  the  parameter  a  on 
going  through  the  critical  value  associated  with  Curve  4,  which  has  a 
bearing  on  the  use  of  Newton's  method  for  solving  Eq.  5b  for  4>  when  the 
coordinates  (a,a^J  are  specified . 

Another  point  worth  noting  is  the  minimum  in  when  a  is 
greater  than  (n2/n^) r  ^  f Curves  5  and  6).  The  minimum  makes  a2(<b) 
stationary,  corresonding  to  a  convergence  of  rays  through  the  point  in 
question .  Such  stationary  points  lie  on  a  caustic  that  extends  from  a  = 
(n2/ni) r i  to  a  -  r,.  By  symmetry  there  are  two  caustics  associated  with 
the  pump  source  at  0^  ,  corresponding  to  positive  and  negative  <p  .  The 
caustics  have  anomalously  high  intensity,  and  are  treated  in  some  detail 
later  in  this  report.  Conventional  ray  tracing  computational  techniques 
that  estimate  intensities  by  counting  the  density  of  rays  tend  to  miss 
the  anomalous  intensities  of  the  caustics . 


Equations  5a  and  5b  are  of  the  form  a  =  af'ifcy  where  a  stands 
for  a1  2  and  a($)  is  the  right  hand  side  of  the  respective  equations. 
Finding  $  for  a  given  point  (a, a)  is  equivalent  to  solving  the  equation 

a(<t> )  -o  =  0.  (7) 


This  then  identifies  the  ray  through  the  point  (a, a)  from  the  source 
located  at  Q„.(or  0.  for  short).  .  It  is  necessary  to  decide  which  of 

Ml  1 

the  two  functional  forms,  5 a  or  5b,  is  the  appropriate  one  to  use,  cor¬ 
responding  to  the  first  or  second  intersection  of  the  ray  with  the 
circle  of  radius  a.  Computationally  this  is  handled  by  generating  a 
table  of  2($)  and  seeing  where  the  target  value  a  lies  in  relation  to 
the  tablulated  values.  If  Newton's  method  is  used,  the  n-th  approxi¬ 
mation  $ generates  the  (n+l)-th  by 


4> 


(n  +  1) 


1  (n)  -  'V 
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where  a  stands  for  the  derivative  da/d$  ,  which  is  given  by 


^  2  2  2  2  2  —  i 

da  /d<t>  =  nR  cos  $  E  (( n .  r.  -  n  R  sin  <b)  2 

1  j=l  J  J~1 


.2  2  2  2  .  2  -  }  . 
(n^  r j  -  n  R  sin  <t>)  } 


da^/dQ  =  da^/d#  -  2ni?  cos# (n2a2  -  n2R2  sin 24>)  2 


The  Pump  Intensity.  (Single  Pass) 


(8a) 

(8b) 


In  the  two-dimensional  approximation  the  intensity  at  the 
field  point  (a, a)  due  to  one  diode  or  one  array  is  estimated  as 
follows.  The  antenna  pattern  at  the  source  is  written  dW/d$  .  (The 
source  is  the  diode  array,  or  its  image.)  At  the  field  point,  consider 
a  line  element  dw  at  right  angles  to  the  ray,  through  which  pass  all  the 
rays  within  the  source  angle  element  d$  .  The  power  radiated  towards 
this  line  is  (dW//d$)d$  and  the  intensity  at  the  field  point  is 
(dW/d$) (d$/dw) ,  multiplied  by  a  factor  that  accounts  for  propagation 
losses.  The  line  element  length  is 

dw  =  | (a2  -  p2 ) 1  da |  (9) 

so  that  the  intensity  is 


I(a, a) 


E  . 
i 


T. 

i 


dW 

d<b(  i) 


da 

d$(i) 


(10) 
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where  all  sources  on  the  circle  are  accounted  for  in  the  summation  over 
i.  The  attenuation  factors  T ^  account  for  reflection  losses  at  the 
boundaries ,  and  absorption  in  the  laser  rod: 


T.  =  exp(-o  s  J  T(6  )  T(Q  )  T(BJ  .  .  .  T(BU  ,)  (11) 

i  -L  f  &  i  4  J  ii*! 


2 

tan  (Q'-Q  ) 

where  T(B1)  =  1  -  - ^ -  (12) 

tan  ( 8'+8j J 


for  P-polarization  of  the  electric  field.  For  S-polarization  tangents 
are  replaced  by  sines,  and  this  is  the  more  usual  case  for  diode  laser 
arrays,  o  is  the  absorption  coefficient  in  the  pump  band  of  the  laser 
rod,  and  s  is  the  length  of  the  ray  from  the  rod  boundary  to  the  field 
point,  given  by 


1,2 


.2  2. 

(a  -  p  ) 


(e\  - 


P2) 


r  2  (a* 


2,i  ,  2 
P  )  (rz 


p2> 1 


(13) 


The  quantities  on  the  right  side  of  Equation  10  are  all 
functions  of  the  respective  variables  $(i)  .  The  angles  8  and  B'  are  all 
expressed  as  inverse  sine  functions  of  a  multiplier  times  sin  $(i)  as 
in  Equation  5a.  Finally,  da /d$(i)  is  given  by  Equation  8.  Computation 
of  the  pumping  intensity  at  the  field  point  thus  comprises  solving  the 
$(i)  in  Equation  7  and  substituting  values  into  Eq.  10. 

Some  Limiting  Cases. 

When  a  equals  p ^  in  the  denominator  of  Equation  10  the  square 
root  factor  becomes  zero  and  the  quantity  a  is  infinity .  However,  it  is 
readily  shown  that  the  denominator  remains  finite: 
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(14) 


|  Lira  {(a2  -  p2)2  dct/d<ti(  i)  }\  =  (n/n^J  R  cos  $(i) 
p  *  a 

At  the  center  of  the  rod  (a=0)  the  denominator  simplifies  to  ( nR/n ,) 
because  $(i)  is  zero.  Both  numerator  and  denominator  simplify  when 
<t>(i)  is  zero  which  occurs  for  a  equal  to  §  or  to  (Q  +  r)  .  In  this  case 
the  transmission  formula  becomes 

m  -  0)  ■  e~','rI  1  al 
-o(r  ¥  a) 


II  { l-(n  .  -  n  .  )2/(n  .  *  n  .  J2} 

j_1  J  J *1  J  J+1 

M-l 

n  t . 
j-i  2 

(15) 


whereas  the  denominator  in  Equation  10  reduces  to 


|  (a2  -  p2)2  da/ d<)>  |  =  |  a  ¥  R(n/n  )  +  a  R  n  {  T.  ( -  ~  )  | 

*  j=l  nj  nj*l  rj 

(16) 

Note  that  Eq.  15  defines  the  quantities  ty  At  a=0  we  get  the  same 
quantity  ( nR/n obtained  earlier  for  the  denominator  and  use  Equation 
15  for  the  numerator .  In  particular  the  intensity  on  the  axis  of  the 
rod  is 


I(axis)  =  (n./nR)  N  T  ( 0,  0) ( dW/d*)  (17) 

1  1  <p=u 

for  N  sources,  where  T±(0,0)  is  the  expression  in  Eq.  15  evaluated  at 
a  =  0  i.e.  the  absorption  factor  is  e  0rl  .  For  a  single  source  (N=l) 
the  last  few  equations  lead  to  a  simple  formula  for  the  intensity  an y- 
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where  along  the  diameter  aligned  with  the  source.  For  example  at 
(r, ,0),  which  is  just  inside  the  laser  rod  on  the  central  ray,  we  have 


I (axis) /!{ r 1> 0) 


(n^r  1/2h) 


(18a) 


where 


M 

(l/2h)  =  Z  (l/n.)(l/r.  -  1/r  .)  .  (18b) 

j  =  2  3  7  3 

It  can  be  shown  that  h  is  the  coefficient  of  the  quadratic  term  in  0^ 
in  the  expression  for  the  optical  path  length  S2  from  the  source  to  the 
edge  of  the  rod 

S2(Q1}  =  S2(0)  *  h9J  +  °(Q1)  (19) 

which  determines  the  curvature  of  the  wavefront  incident  on  the  laser 
rod  surface. 

Normalization  of  the  Pump  Intensity. 

In  the  two-dimensional  geometry  the  antenna  pattern  from  the 
pump  source  is  considered  to  be  a  function  of  the  single  variable  $  . 
The  quantity  dW/dg  in  Eq.  10  needs  to  be  expressed  in  terms  of  P/L,  the 
power  per  unit  length  radiated  by  the  diode  array  of  length  L.  It  is 
clear  that 

P/L  =  J  (dW/d$)  (20) 

and  it  is  instructive  to  assume  that  the  pattern  dW/d $  is  a  Gaussian 
function: 
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dW/dip  =  (dW/d$) 


(21) 


.  A  exp  - r 4»/4>  )' 
g=u  e 


Then  the  integral  in  Eq.  20  is  (dW/db)  A  it  2  which  evaluates 

$=0  e 

(dW/d$).  .  so  that  Eq.  21  becomes 
0  =  0 


dW/db  =  (P/L)  (1/b  it'  J  exp  -  r )‘ 
e  e 


Thus  it  is  easy  to  substitute  for  dW/d$  in  Eq.  10,  or  for  (dW/dQ)  in 

<J>=0 

Bg.  17  to  calculate  the  intensity . 


In  some  of  our  computational  illustrations  the  diode  laser 

2 

antenna  pattern  has  been  normalized  simply  to  exp  -  ($/<()  )  .  When  this 

e  ± 

is  done  the  power  radiated  per  unit  length  of  diode  array  is  $  i r1  ,  and 

e 

the  computed  intensities  in  the  rod  refer  to  this  source  strength.  For 

unit  source  strength  (1  watt /cm)  the  computed  intensities  would  need  to 

_1 

be  divided  by  $  v*  .  These  cases  will  be  identified  in  the  figure 
e 

captions  where  they  occur.  It  is  more  convenient  to  build  the 

s 

normalization  factor  l/('\>e it'  )  right  into  the  computer  program  as  we 
have  done  more  recently,  so  that  the  results  give  intensity  (watts /cm2 ) 
corresponding  to  unit  source  strength  (1  watt / cm) .  With  N  source 
arrays,  the  pump  light  near  the  axis  is  arrriving  in  N  different 
directions .  The  units  of  intensity  that  we  use  ( watts /cm^J  are  still 
valid  because  the  total  transition  probably  for  pumping  an  atom  or  ion 
is  the  sum  over  the  individual  contributions  irrespective  of  their 
directions . 


A  Computational  Example. 


The  following  example  is  used  to  illustrate  the  use  of  the 
theory  so  far  presented.  Consider  a  YAG  rod  of  radius  0.3175  cm.  (This 
is  a  quarter  inch  diameter.)  The  surface  of  the  rod  is  assumed  to  be 
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smooth  so  that  the  laws  of  refraction  apply.  The  rod  is  surrounded  by 
cooling  water  enclosed  in  a  glass  jacket.  Five  pumping  sources  are 
equally  spaced  on  a  circle  of  radius  4  cm.  We  are  interested  in  the 
pump  intensity  distribution  inside  the  rod,  which  is  assumed  to  have  an 
attenuation  factor  of  3,  so  that  the  total  attenuation  along  a  diameter 
is  exp  -  (0.3175  x  3)  =  0.39  .  Light  that  passes  through  the  rod  is 
wasted  in  this  computational  example.  The  pumping  sources  are  assumed 
to  have  equal  to  6.46  degrees.  (The  source  strength  is  0.1998,  which 

i  e 

is  g  irJ  .  )  These  sources  would  most  likely  be  the  virtual  images  of 

e 

real  diode  arrays  situated  on  a  smaller  radius  than  4  cm,  radiating  into 
a  positive  lens,  because  the  arrays  normally  have  a  wider  antenna  pat¬ 
tern.  The  pump  light  is  assumed  to  have  S-polarization  of  its  electric 
field  i.e.  polarized  along  the  direction  of  the  array.  The  fraction  of 

pump  light  radiating  at  angles  less  than  <(>  is  given  by  the  Error 

max 

Function  erf  (0.940).  (This  number  is  0.816.)  Light  at  larger  angles 
is  wasted,  as  it  cannot  enter  the  laser  rod. 

At  any  fixed  radius  a  the  intensity  has  five  equal  cycles  in 
the  360  degrees,  and  shows  symmetry  about  any  diameter  through  a 
source.  Therefore  we  need  to  tabulate  only  within  the  range  0  to  36 
degrees . 


In  the  computations  the  refractive  indices  assumed  are  1.823, 
1.333,  1.5  and  1  for  the  media.  The  geometry  is  shown  in  Figure  5. 

Figure  6  shows  intensity  versus  angle  at  radial  steps  of 

r^/10.  Notice  that  from  the  center  (a/r^  =  0)  out  to  a  =  0.7r ^ 

there  is  very  little  angular  variation  in  intensity .  At  radii  exceeding 

p  (0.7312  r  )  the  effects  of  caustics  show  up.  ( See  lowest  three 

m  1 

curves  in  Fig.  6.)  There  is  a  smooth  drop  of  intensity  with  increasing 
radius  along  radii  bisecting  or  passing  through  the  sources  because  in 
these  special  cases  no  caustics  cross  these  radii. 
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ALPHA 

C88-8952 

Figure  6.  Intensity  versus  angle  at  various  radii  (normalized  to  rod 
radius  r^)  for  the  geometry  shown  in  Figure  5.  The  antenna 
width  $  is  6.46  degrees  (0.1127  radian),  so  that  the  power 
radiated  from  each  source  is  0.1998  watt  per  centimeter. 


General  Considerations . 


Irrespective  of  the  distance  R  (locating  the  point  source) 
there  are  two  segments  unilluminated  as  illustrated  in  Figure  7.  These 
segments  are  bounded  by  arcs  of  angle  (a^  -  a  )  evaluated  at  a  =  r^, 
i .  e. 


Arc  =  it  -  2  sin  (n^/n^J 


Likewise  the  perpendicular  distance 


Pm  =  n2ri/ni 


is  invariant .  What  varies  with  R  is  the  angle  x  that  p  makes  with 

m  *m 

the  central  ray.  This  angle  is  given  by  the  following  equation. 


’  *  ,ai  '  V«u 


M 


.  -1 


=  0  +  it  +  Z  sin  (n.r./n  ,r  .  ,) 
j-3  2  1  1  I'1 

M 

-  fml  sin1(n2rl/n.r.) 


(23) 


The  quantity  nR  should  be  chosen  to  set  at  the  most  desirable  value 
with  respect  to  the  antenna  pattern  of  the  diode  array.  If  the  angular 
width  of  the  pattern  is  excessive  a  lens  might  be  used  to  alter  it. 
Another  consideration  is  the  variation  in  pump  intensity  due  to  the 
focussing  property  of  the  optical  arrangement.  In  terms  of  x  measured 
along  the  central  ray  as  illustrated  in  Figure  7,  the  intensity  obtained 
from  Equations  10,  15  and  16  is 
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I(x) 


1  “l.  ~01i 

=  R  ^  e 


J¥-l 

n 

j=3 


t  .  jfx; 


where  the  t- functions  are  those  in  braces  in  Equation  15,  and  J(x)  is 
given  by 


J(x)  =  e°  X/| (1  +  hx> 


where 


(24b) 


with  h  defined  in  Eq.  18b.  (Equation  24  applies  on  the  x-axis  only.)  b 
may  be  made  positive  or  negative  by  appropriate  choice  of  parameters. 
It  is  readily  shown  that  dJ/dx  is  zero  and  d2J/dx 2  is  minimum  when  b 
equals  a  ,  which  thus  gives  the  condition  for  the  most  uniform  pump 
distribution  along  the  central  ray  near  the  center  of  the  rod,  from  a 
single  array.  (This  is  not  a  really  significant  result,  though.) 

The  following  is  an  approach  to  the  problem  of  seeking  the 
most  uniform  pump  distribution  around  the  axis  of  the  rod  for  multiple 
diode  arrays.  The  transverse  variation  is  best  approximated  near  the 
center  by  looking  at  the  values  along  the  Locus  p  =  a,  where  from 
Equations  10  and  14  we  have  for  a  single  array. 


I( $  =  sin 


V 

nR 


n 

j 


«v  35 


—  R  cos b 
n  1 


(25) 
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wherein  p  will  be  replaced  by  x.  (Here  we  assume  that  the  source  lies 
on  the  y -axis,  so  that  p  lies  approximately  along  x  for  small  p.)  This 
expression  is  an  even  function  of  x  that  will  be  expanded  as  far  as 
quadratic  terms.  The  various  factors  give: 


2  2  i 

■ o(r '  -  p  )  1 


e~°r 1 (i  +  a  x2} 


dW  , . _  _  .  .  1  .22 

—  ( if  Gaussian)  =  1  -  (  — - — )  x 
d< J>  n  R$ 


_J _ ,  .  1  .2  2 

cos<t>  *  2  n  R  x 


Neglecting  the  variation  in  the  T(B)  functions,  the  result  is 


I(x)/I(o)  -  1  *  x2{  /  f-  -  )2( ♦  0(x)4 

r  7 

1  ^  e 


,  .  o  .  1  .2.  2 

1  *  (  i  7 - f  p  7  )  >  x 

r,  n  R  <p 
I  e 


(26) 


where  we  assume  <<  2. 

e 

Consider  four  equally  spaced  pump  beams  (two  opposing  pairs  at 
right  angles) .  Along  one  of  the  axes  the  variation  from  one  pair  of 
beams  goes  like 
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\j(x) I  +  ]jf-xj|  =  2(cosh  ax  -  bx  sinh  ax) /(l  -  b 2  x2 ) 


if  bx  is  less  than  unity,  while  the  contribution  of  the  other  two  is 
twice  the  right  side  of  Equation  26  above.  The  coefficient  of  x 2 
vanishes ,  giving  the  condition  for  most  uniform  intensity  variation 
around  the  axis,  when 


where  b  =  b(R)  as  in  Equations  (24b)  and  (18b).  Once  the  laser  rod  and 
water  jacket  sizes  are  chosen  as  in  Figure  2  for  example,  Equation  28 
can  be  represented  by  a  set  of  curves  in  the  ($e,R)  plane  with  parameter 
a  .  Figure  8  shows  such  a  set  for  a  YAG  rod,  water  and  glass  water 
jacket  of  dimensions  indicated  in  the  figure. 


Although  the  conditions  for  uniform  pumping  have  been  derived 
for  four  beams,  obviously  the  same  conditions  would  apply  for  eight  or 
twelve.  We  can  assume  that  they  would  be  accurate  for  any  other  number 
above  four. 
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88-7503 


Figure  8.  Antenna  width  in  formula  exp  [-($/$  )2]  that  gives  flattest 
pump  distribution  of  axis  of  rod  versus  diode  radial 
distance  for  various  assumed  pump  band  coefficients  o  in 
the  expression  exp  [-ox]  .  The  geometry  in  shown  in 
Figure  5. 
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The  motivation  for  seeking  a  uniform  energy  deposition  over 
the  cross  section  of  the  rod  is  that  it  comprises  a  first  step  in  under¬ 
standing  the  more  complete  problem  that  optimizes  distribution  of  ther¬ 
mally  induced  strain,  population  inversion,  and  temperatures . 

Multipass  Analysis . 

The  results  shown  in  Figure  6  show  that  the  pumping  intensity 
can  be  made  fairly  uniform  over  about  one  half  the  cross  section  area, 
but  there  is  a  pumping  efficiency  penalty  in  that  a  large  fraction  of 
the  pump  light  is  transmitted  through  the  rod.  The  efficiency  can  be 
increased  by  increasing  the  absorption  coefficient  o  .  Another  approach 
is  to  reflect  the  transmitted  light. 

Suppose  that  the  pump  arrays  are  separated  by  sections  of  a 
cylindrical  mirror  of  radius  R.  Then  a  ray  at  angle  $  ,  after  passing 
through  the  laser  rod  strikes  the  mirror  at  angle  of  incidence  equal  to 
<(>  ,  because  as  shown  in  Figure  9,  the  ray  path  is  symmetrical  about  the 
perpendicular  p.  After  reflection  the  ray  will  leave  the  mirror  surface 
at  an  angle  $  .  All  rays  that  pass  through  the  rod  will  make  a  set  of 
rays  displaced  an  angle  2 t  from  one  set  to  the  next,  where 

r  =  i  (ct1  +  <*2)  (23) 

(t  is  negative  when  is  positive .  See  Figure  9.)  The  N  arrays  lay 

down  a  pattern  of  2N  rays  belonging  to  angles  ±  4>  on  the  first  pass. 

These  rays  lie  within  the  ring  bounded  by  radii  p(b)  and  r .  After  one 

reflection  of  each,  there  are  now  an  additional  2N  rays  within  the  same 

ring;  and  each  successive  reflection  produces  a  further  2N  rays .  Light 

at  a  smaller  angle  similarly  produces  a  pattern  of  rays  bounded  by 

pO'j  and  r  ,  overlapping  the  first  ring  that  we  considered .  The 

illumination  within  the  circle  a($  <  $  )  has  been  found  by  computation 

n 
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Figure  9.  In  one  pass  through  the  system  a  ray  from  the  source  at  angl 
<t>  advances  in  position  on  the  circle  of  radius  R  by  2t  . 

A  plot  of  2t  against  <t>  for  a  typical  case  is  shown  in 
Figure  11. 


to  be  essentially  independent  of  angle  a  of  the  (a, a)  ,  for  five  or  more 
arrays.  The  addition  of  interlaced  rays  in  the  same  domain  produced  by 
the  reflections  from  the  cylindrical  mirror  would  not  be  expected  to 
essentially  change  this  circular  symmetry .  We  assume  then,  without 
mathematical  proof,  that  the  intensity  in  the  rod  will  be  at  least  as 
symmetrical  as  for  the  single  pass  and  that  the  contribution  from  all 
rays  at  any  angle  0  will  be  equivalent  to  the  single  pass  contribution 
from  an  antenna  intensity  dW' /d$  related  to  dW/d$  by 


dW'  dW  _  .  2 ...  3 

^  ^  [i  +  P($)  +  P  ($)  +  P  (t)  *  •  •  •  •] 


ft'U  -  Pit)  i 


(30) 


where 


-2o[r‘  - 


2  2 

n  R  .  2A1  / 
— —  sm  iti] 


Q(b)  * 


II  T  ( B  J  e 

j  J 


(31) 


Here  the  T(B)  are  the  transmission  factors  given  in  Equation  12,  the 
exponential  factor  accounts  for  single  pass  absorption  in  the  rod,  and  G 
is  the  reflection  coefficient  of  gold  (0.974).  The  radial  distribution 
for  the  multipass  case  is  thus  estimated  by  modifying  the  antenna  pat¬ 
tern  as  in  Equation  30,  then  doing  the  single  pass  computation. 


Using  Equations  30  and  31  we  have  calculated  the  radial 
distributions  shown  in  Figure  10,  where  the  diode  arrays  are  at  R  =  4 
cm,  and  dW/d$  is  chosen  to  be  Gaussian  with  equal  to  6.7  degress. 
With  the  absorption  coefficient  o  equal  to  3  the  parameters  are  set  to 
the  point  X  in  Figure  8  which  predicts  a  uniform  pumping  near  the 
axis.  The  calculated  distribution  in  Figure  10  suports  this  prediction 
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a  (cm) 


88-7504 


Figure  10.  Radial  distributions  for  near  optimum  parameters  aimed  at 
uniform  distribution.  (Multipass  calculation.)  These 
distributions  are  computed  for  P-polarization,  and  are  not 
normalized  to  unit  power  per  unit  length  from  the  source. 
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fairly  well,  especially  in  view  of  the  fact  that  the  reflective  losses 
at  the  glass  and  YAG  surfaces  were  not  folded  into  the  angular  distri¬ 
bution  of  the  pump  in  the  optimization  analysis .  (The  curves  in  Figure 
10  were  calculated  along  the  line  a  =  0  i.e.  along  a  radius  aligned 
with  one  of  the  seven  arrays.  If  the  calculations  had  been  performed  at 
any  different  a  ,  the  differences  would  have  been  undetectable  on  the 
curves  until  a  reached  0.24.  The  observational  results  would  probably 
show  more  symmetry  than  the  calculations.) 

The  rays  that  leave  an  array  spread  over  a  significant  arc  on 
arriving  at  the  reflecting  surface,  as  illustrated  in  Figure  9.  The 
spread  is  shown  in  Figure  11  for  a  typical  geometry .  If  is  small, 

say  2  degrees,  the  energy  transmitted  through  the  rod  would  be  concen¬ 
trated  over  a  range  of  about  25  degrees.  The  dependence  of  this  range 
on  geometry  has  not  been  investigated .  If  the  range  is  small  it  makes 
sense  to  have  an  odd  number  of  arrays  if  they  are  arranged  inside  the 
circle  R  (and  imaged  at  R) .  Alternatively,  the  diode  arrays  might  be 
outside  the  circle  R  and  imaged  onto  slits  in  the  reflector.  The  choice 
between  mounting  the  diode  inside  the  reflector  (with  lenses  forming 
virtual  images)  and  outside  (with  real  images  at  the  slits)  would  prob¬ 
ably  be  made  on  the  basis  of  mechanical  convenience  and  simplicity 
(favoring  outside  mounting) . 

The  conditions  for  uniform  pumping  in  the  region  of  the  axis 

of  the  rod  as  expressed  in  Eq.  28  and  Fig.  8  is  modified  somewhat  in  the 

multipass  case,  owing  principally  to  the  effective  broadening  of  the 

antenna  pattern  in  Eq.  30.  Again  a  Gaussian  antenna  pattern  exp 
2 

-($/$  )  will  be  used  to  illustrate  the  case,  and  again  the  relevant 
e 

expressions  in  $  and  x  will  be  expanded  as  far  as  second  order  in  these 
variables  with  x  equal  to  (nR/n^)<b  .  Expansion  of  the  right  side  of  Eq. 
31  gives 
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Figure  11 


.  Angular  positions  of  rays  after  one  pass.  The  geometry  is 
that  shown  in  Figure  5. 


p  =  p  +  pi 
O  2 


r2  -2"l 
p  =  GT  e 
o  o 


2  2  ~ 2ar 1 

and  p_  =  G{2T  T  +  T  (nR/n  )  (c/r  )}  e 

2  O  2  O  1  1 


In  these  expressions  Tq  and  T^  come  from  the  expansion : 


n  T(  B  .)  =  T  +  TA 
i  1  °  2 


where  T  =  II  t  .  , 
°  1=1  3 


9  ?  2 

T„  =  -2  n  R  T  n  (1  -  t.)/(t.n.n.  .  r  .) 
2  O  .  7777  +1  ] 

J=1 


(See  Eg.  15  for  the  definition  of  tj-) 


This  expression  for  T2  is  valid  for  S-polarization.  It  must  be  multi¬ 
plied  by  the  factor  (-2)  for  P-polarization. 

The  antenna  pattern  in  Eg.  30  to  second  order  in  $  is: 


dW 

d<t> 

1-0. 


—)  t2} 


1  -  0 


TTC-1529-R 


30 


Considering  four  pump  beams  as  in  Eqs .  27  and  28,  we  find  the  expansion: 


I(x) /I(o) 


1 


+ 


4 

~  [(a  -  b) 


2 


+  b 


2n ‘ 


2  2 
n  R 


-  i  )) 


(35) 


Thus  the  left  side  of  Eg.  28,  which  gives  the  condition  for  most  uniform 
pump  intensity  variation  near  the  axis,  is  to  be  replaced  by 


(36) 


Here  the  last  term,  }  ,  comes  from  cosij>  in  the  denominator  of  Eg.  25, 

T2/'To  comes  from  the  ( H-l )  factors  T($)  in  the  numerator,  and  the  term 

P^/n  “  PQ)  comes  from  the  multipass  effect  in  Eg.  34.  The  results 

summarized  in  Figure  8  for  example  can  therefore  be  readily  extended  to 

the  multipass  case  by  reading  the  vertical  axis  as  ,  and  subsequently 

e 

evaluating  4>e  by  substituting  the  values  of  pQ,p^, ro  and  T2  into  Eq. 
36. 


The  optimum  pumping  intensity  distribution  in  terms  of  mini¬ 
mizing  thermal  equilibration  times,  or  amplifying  a  single  (Gaussian) 
mode  input  signal,  would  correspond  to  the  coefficient  of  x*  in  Eq.  35 
being  negative .  The  corresponding  relation  amongst  the  design  para¬ 
meters  may  be  readily  found  and  expressed  as  a  set  of  curves  like  those 
of  Fig.  8. 
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Application  of  the  Theory. 


One  way  to  use  the  theory  above  is  to  consider  a  laser  rod  and 

water  jacket  of  specific  radii,  such  as  those  values  shown  in  Figure  5; 

and  then  to  regard  R  and  4>  as  adjustable  parameters  in  an  exercise 

e 

aimed  at  producing  a  relatively  uniform  pumping  intensity  profile  in  the 
middle  of  the  rod.  This  implies  that  whatever  the  antenna  character¬ 
istic  width  is  in  the  diode  laser  arrays  chosen  for  pumping,  one  can 
modify  it  to  the  desired  value  by  use  of  a  lens,  and  place  to  at  the 
desired  radius  by  straightforward  design.  Figure  12  illustrates  poss¬ 
ible  methods  for  this  modification.  Eq.  28  (with  the  left  side  replaced 
> 

by  <)>e  as  defined  in  Eg.  36)  may  be  regarded  as  a  functional  relation 

between  R  and  the  optimum  value  of  $  for  various  values  of  o  .  The 

e 

quantity  b  captures  the  focussing  effects  of  the  media  from  the  pump 
source  to  the  rod.  The  equations  have  been  programmed  in  BASIC  on  our 
IBM  PC.  OPTIM.SIN  and  OPTIM.TAN  give  the  optimum  4>e  for  given  values 
of  R  for  S  and  P  polarizations  respectively.  An  example  of  the  results 
is  given  in  Table  I  for  the  geometry  of  Figure  5  (with  variable  R) , 
using  S-polarization.  (Fig.  8  contains  the  same  kind  of  data.)  The 
reflectivity  of  gold  is  about  0.974  at  0.8  micron  wavelength.  The 
computations  have  assumed  a  wall  reflectivity  of  0.95,  i.e.  roughly 
twice  the  loss  of  a  perfect  gold  coating  on  a  polished  surface.  The 
effective  reflectivity  is  reduced  by  the  presence  of  the  entry  slits  if 
external  focussing  optics  are  used,  or  by  the  internal  optics.  If 
internal  optics  are  used  in  conjunction  with  reflecting  walls,  these 
optical  components  should  obviously  be  very  small. 

The  angles  tablulated  as  4 >^(3  or  5)  are  the  characteristic 
antenna  widths  of  optimum  uniformity  (according  to  Eqs.  28  and  36)  near 
the  axis  of  the  rod  in  cases  where  a  is  3  or  5  respectively .  It  is 
seen  that  higher  absorption  coefficient,  which  tends  to  lower  the 
illumination  on  the  axis  compared  to  that  at  the  surface,  requires 
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Figure  12.  Adjustment  of  source  position  radius  R  and  antenna  width  $ 
by  use  of  internal  or  external  optics.  e 


TABLE  I 


OPTIMIZATION  OF  PARAMETERS 


The 

Geometrical  Parameter  b. 

Optimum  $  , 

e 

and  Error  Function  Argument 

foe 

Various  Values 

of  R.  The 

Fraction  of 

Power  Radiated  in 

Angles  less 

than 

the  Critical  Angle  is  given  by  Erf(x) 

The  Absorption 

Coefficient 

is  3  cm~ 1  in  x(3), 

of  the  Wall  is  0.9 5 

<J>  (3)  and 
e 

5  cm  1  in  x( 5),  <p  (5)  .  The 

e 

Reflectivity 

R 

b 

*J3) 

e 

x(3) 

5) 

e 

x(5 ) 

1 

0.0349 

31.97 

.  783 

22.68 

1.104 

4 

-1.332 

6.  74 

.901 

4.91 

1.236 

10 

-1.606 

2.57 

.  944 

1.90 

1 . 279 

50 

-1.  752 

0.  500 

.  970 

.372 

1.304 

5  00 

-1.  784 

0.0497 

.  975 

.0370 

1 . 310 

-50 

-1.825 

0.493 

.  982 

.  368 

1.317 

-10 

-1.970 

2.40 

1.01 

1.81 

1.343 
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as  expected .  The  pumping  efficiency ,  i.e.  the 


smaller  antenna  width  <p 

e 

fraction  of  diode  laser  power  absorbed  in  the  rod,  contains  a  factor 
that  accounts  for  the  fraction  of  radiation  at  angles  less  than  the 
critical  angle,  which  is  the  error  function  of  <t>  (critical )  /  <s>  .  This 

ratio  is  denoted  by  x  in  the  tables.  Both  the  critical  $  and  <p 

e 

decrease  with  increasing  R (effective  source  radius)  but  x  increases . 
From  the  tablulated  values  for  R  =  5 00  cm  it  is  seen  that  a  collinated 
beam  is  better  than  a  diverging  one;  and  from  the  values  for  R  =  -50  and 
-10,  it  is  clear  that  a  converging  beam  is  better  still  in  terms  of  this 
efficiency  factor.  However,  the  internal  optics  make  assembly  and 
adjustment  more  difficult. 

The  approach  that  we  have  used  for  deducing  the  conditions  for 

most  uniform  pumping  near  the  axis  is  based  on  a  quadratic  approximation 

for  the  intensity .  Obviously  the  assumed  Gaussian  antenna  pattern  has  a 

more  rapid  fall-off,  and  the  result  is  that  the  parameters  deduced  from 

the  theory  tend  to  give  a  maximum  on  axis  as  shown  in  Fig.  6  and  the 

upper  curve  in  Fig.  10.  If  a  very  uniform  pumping  intensity  is  desired 

the  approach  to  be  used  is  that  suggested  by  Fig.  10,  which  is  to  s  rect 

a  value  of  <)>  that  belongs  to  a  somewhat  smaller  a  (absorption  co- 
e 

efficient) .  However  this  lowers  the  value  of  x.  The  value  of  Erf ( 1 . 16) 
is  0.90,  so  that  x  should  be  greater  than  1.16  if  this  efficiency  factor 
is  to  be  at  least  90  percent.  (As  in  most  laser  applications,  while 
mathematical  calculations  are  simplified  by  the  assumption  of  a  Gaussian 
intensity  profile,  a  more  desirable  profile  is  one  that  has  a  flatter 
top  and  steeper  sides.) 

Computer  programs  for  determining  the  pump  distribution  are 
readily  developed .  The  theory  is  simple  enough  for  programs  to  be 
written  on  a  hand-held  calculator ,  such  as  the  Hewlett  Packard  HP 
41CV.  The  calculator  is  slow  but  it  is  has  the  advantage  of  " double 
precision"  built  in  automatically,  which  increases  the  reliability  of 
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the  programs  against  crashes.  Our  IBM  PC  BASIC  program  "AUTOPUMP" 
generates  the  pump  intensity  distribution  within  the  basic  sector  at 
eleven  radii  (including  zero),  and  eleven  angles  at  larger  radii  and  six 
angles  at  the  lower  radii  where  the  angular  variation  is  smaller.  The 
program  can  handle  up  to  10  concentric  media  (and  more,  as  modified  by 
Pharos).  This  program  and  its  variants  have  been  exercised  extensively 
and  many  results  that  should  be  obvious  a  priori  are  revealed .  Some  of 
these  observations  are  as  follows. 

The  inner  boundary  of  the  caustics  is  at  radius  a  =  n  r  /n  . 

2  2  11 

The  area  inside  this  radius  contains  ° ^  cross  section,  and 
one  can  make  the  intensity  of  pump  light  quite  uniform  in  this  region. 
Some  typical  values  of  the  cross  section  area  ratios  are  given  in  Table 
II.  The  relative  non-uniformity  in  the  outer  part  of  the  rod  has  two 
principal  causes,  the  caustic  and  the  proximity  of  the  sources.  As  an 
extreme  example  of  the  latter,  if  the  sources  are  close  to  the  rod  and 
have  narrow  antenna  patterns,  and  if  in  addition  the  absorption  through 
the  rod  is  large,  then  the  pump  intensity  shows  a  set  of  strong  peaks 
around  the  outer  portion  of  the  rod  each  facing  a  pump  source. 
Obviously,  increasing  the  number  of  sources  on  the  circle  improves  the 
uniformity .  At  least  three  sources  would  seem  to  be  necessary .  Four  or 
five  would  be  better. 

Pumping  Efficiency . 

The  light  emitted  by  the  diode  lasers  is  absorbed  in  the  laser 
rod,  the  reflecting  walls  of  the  laser  head,  and  in  the  cooling  fluid 
(e.g.  water)  and  its  jacket  (e.g.  glass  of  fused  quartz).  The  absorp¬ 
tion  in  pure  water  is  about  0.023  cm~^  at  0.800  micron  wavelength,  and 
will  account  for  about  one  percent  of  the  pump  lignt  absorption .  Of  the 
total  pump  light  available  a  fraction  n  is  absorbed  in  the  laser 
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TABLE  II 


Medium  I 

(n^ 

YAG 

(1.823) 

YAG 

YAG 

YAG 

YLF 

(1.634q, 

YLF 


Area  with  Smooth  Pump  Intensity  -  Inside  the 


Medium  II 


(n2) 


Water 

(1.333) 

Ethylene  Glycol 
(1.427) 

Quartz 

(1.453) 

50-50  Eth.  Gly. -Water 
(1.380) 

Water 

1.631e) 


50-50  Eth.  Gly. -Water 


Caustics 


Area  Ratio 
(n2/nL)2 


0.  53 


0.61 


0.66 


0.57 


0.67 


0.  71 


TTC-1529-R 


37 


rod.  The  quantity  n  is  called  the  pumping  efficiency  here,  and  it  can 
be  calculated  with  accuracy  consistent  with  that  of  the  pumping  inten¬ 
sity  calculations,  i.e.  limited  by  the  assumptions  that  go  into  Eqs.  30 
and  31  in  the  multipass  analysis .  There  are  two  simple  approaches  to 
the  calculation  of  n  .  One  is  to  figure  out  the  absorptions  of  a  single 
ray  as  it  makes  successive  passages  through  the  rod.  The  other  is  to 
integrate  al  over  the  volume  of  the  rod. 

The  successive  absorptions  in  the  rod  and  in  the  surface  of 

the  reflector  at  radius  R  respectively  can  be  added  using  the  multipass 

analysis .  The  notation  used  in  Eq .  31  is  simplified  by  expressing  the 

-aS(4>) 

exponential  factor  as  e  and  by 

T(<t>)  =  T(  BI;  T(&2)  T(&3)  .  .  . 

The  absorbed  powers  are: 

r  a  -1  -nS 

P(rod)  =  J  (1-p)  T(l-e  )(dW/d*)d *  (37) 

-4 

Ja  -7  y  _  rj  S 

(1-p)  T*(l-G)  e  ( dW/d$)d$  (38) 

-4 

These  formulae  count  only  the  light  radiated  at  angles  less  than  4  ,  the 

antenna  pattern  cut-off  or  ,  whichever  is  smaller.  If  <(>  is  smaller, 

the  pumping  efficiency  contains  the  factor  W($  ) /W(Q)  (assuming  that  W 

M 

is  symmetric)  which  accounts  for  pump  light  that  cannot  reach  the  laser 
rod.  Pump  light  reflected  or  scattered  at  the  interfaces  between  suc¬ 
cessive  media  has  been  neglected  in  Eqs.  37  and  38,  thus  representing  a 
deficit  of 

P (neglected)  =  J®  (l-p)"1 ( 1-T) ( l+Te~°S ) (dW/d$)d$  (39) 

-4 
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Some  of  this  will  find  its  way  into  the  rod. 


From  Eg.  34  it  is  clear  that  the  pump  light  intensity  at  the 

center  of  the  rod  is  increased  by  the  factor  (1  -  p )  ^  by  providing  a 

reflecting  surface  at  radius  R.  The  total  increase  in  pumping  due  to 

the  reflecting  surface  is  found  by  computing  Eq.  37  with  the  appropriate 

functional  form  for  p(b)  and  dividing  by  the  same  integral  evaluated 

with  p  equal  to  zero:  The  result  is  somewhat  greater  than  (1  -  p  )  1 

o 

because  the  function  dW' /d$  is  broader  than  dW/d $  as  seen  in  Eq.  34; 
but  the  factor  (1  -  pq)  *  is  a  good  approximation.  (Recall  that  pQ  is 
the  total  attenuation  along  the  central  ray.) 

Consider  the  following  example: 

n i  =  1.923,  r =  0.3175  cm,  a  =  3  cm  1  (YAG) 

n2  =  1.333,  r2  -  0.5  (Water) 

n^  =  1.5,  r3  ~  0.65  (Glass  water-jacket) 

n  =  1,  R  =  4  cm  (Location  of  diodes) 

dW/d$  =  exp  -(ty/6.7  deg.)2 

<b  =  6.073  deg. 

PI 

G  =  0.974  (reflectivity  of  gold) 

Pump  Polarization:  P 

T(o)  =  .9336  (Air-glass-water-YAG) 
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1.1446  . 


Then  p  =  0.1263  ,  and  ( 1  -  p  ) 
o  o 

Of  the  diode  power  radiated  into  angles  less  than  6.073  deg.  the  frac¬ 
tion  0.9365  is  absorbed  in  the  YAG  rod.  If  G  *  0  (single-pass  pumping) 
the  fraction  absorbed  is  0.7773.  The  pumping  is  increased  by  the  factor 
0.9365/0.7773  which  is  1.2048,  which  is  to  be  compared  to  the  simple 
approximation  that  gives  1.1446.  Notice  that  exp( -2or  )  is  0.1488  while 
the  corresponding  quantity  for  the  critical  ray  is  0.2727.  The  pumping 
efficiency  0.9365  is  now  multiplied  by  W(b  )/W($)  which  is  the  error 
function  erf (6.073/6.7)  =  0.800.  Multiplying  this  by  0.9365  drops  the 
efficiency  to  0.  75.  .  On  the  other  hand,  much  of  the  drop  from  unity  to 
0.9365  is  due  to  having  T(o)  =  0.9336.  (T  increases  with  increasing 
4)  because  of  the  polarization  chosen.)  Without  Fresnel  losses  the 
fraction  absorbed  would  go  from  0.9365  to  0.9945.  Thus  the  Fresnel 
reflection  losses  account  for  5.8  percent. 

In  summary,  the  various  pumping  efficiency  factors  in  the 
chosen  example  are  as  follows. 

Radiated  into  4>  <  4>„  (Gaussian) :  0.800 

M 

Fresnel  for  P-polarization :  0.9417 

Absorption  in  rod  (multi-pass) :  0.9945 

Net  efficiency :  0.75 

Multi-pass /Single  pass  1.20 

The  net  efficiency  should  be  multiplied  by  the  approximate  additional 
factors  0.99  to  account  for  absorption  in  the  cooling  water,  and  0.95  if 
the  polarization  is  S  instead  of  P.  These  will  bring  the  net  efficiency 
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down  to  0.70.  The  biggest  loss  factor  is  that  due  to  the  lost  power  in 
the  Gaussian  wings  of  the  antenna  pattern  of  the  diodes,  and  many  other 
numerical  examples  have  verified  this  same  effect. 

The  factor  G  that  occurs  in  Eg.  38  and  in  Eq.  31  for  p  (which 
occurs  in  Eqs .  37  and  38)  refers  to  the  reflectivity  of  the  wall,  but 
also  accounts  for  the  obscuration  of  the  reflecting  wall  by  any  internal 
optics.  As  an  example,  if  we  choose  R  =  5 00  cm  this  implies  that  ve 
have  internal  optics  in  front  of  the  diodes  as  shown  in  Fig.  12.  A 

reasonable  value  to  use  for  G  in  this  case  would  be  0 . 9.  Since  the 

reflectivity  of  clean  gold  is  0.974  at  0.8  micron,  a  value  of  G  higher 
than  0.9  can  be  assumed  for  R  =  4  cm,  especially  if  the  diode  lasers  are 
outside  the  gold  reflector,  with  the  pump  light  entering  through 
slits.  This  case  also  is  shown  in  Fig.  12.  Because  the  effective 

value  of  G  is  larger,  R  =  4  cm  gives  a  higher  n  than  R  =  5 00  cm.  The 

intensity  plots  of  the  R  =  50  cm  case  are  shown  in  Figure  13. 

We  have  assumed  in  these  computations  that  the  light  reflected 
from  the  interfaces  is  lost,  whereas  it  is  reflected  back  from  the  walls 
towards  the  rod  again.  Eq.  31  is  readily  modified  to  account  for  this 
additional  feedback.  The  result  is: 

P (*)  =  (1  -  t  *  te~°S)G  (40) 

where 

t  =  n  T(B  .) 

J  J 

2  2 

,  „  2  n  R  .  2  ,  ,  { 

and  S  =  2[r ^  -  — —  sm  4>] 

"l 

In  order  of  importance,  we  can  list  the  measures  recommended 
for  high  pumping  efficiency . 
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Figure  13.  Pump  Intensity  distributions  for  R  =  50  cm  at  different 

angular  positions,  corresponding  to  the  "Virtual  Image"  case 
in  Fig.  12.  Equally  uniform  distributions  can  be  generated 
by  computation  of  the  other  two  cases  with  optimization  of  <t>. 
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(i)  Use  a.  highly  reflecting  wall  for  pumplight  feedback. 


(ii)  Use  optics  external  to  the  reflecting  wall  for  control¬ 
ling  the  pump  beam  characteristics ,  and  have  the  pump  light  enter 
through  narrow  slits  in  the  wall. 

(Hi)  Design  the  optics  so  as  to  bring  light  from  the  wings 
of  the  antenna  pattern  into  the  central  portion,  so  that  it  can  enter 
the  rod. 


(iv)  Apply  anti-reflection  processing  to  water-jacket, 
lenses,  and  other  surfaces  intersecting  the  pump  rays. 

(v)  Minimize  the  number  of  concentric  media  in  the  design  so 
as  to  reduce  losses  from  the  pump  to  the  laser  rod. 

(vi)  As  an  experimental  matter,  check  the  absorption 
coefficient  of  the  pump  light  in  the  cooling  liquid. 

The  second  method  of  calculating  the  efficiency  is  outlined 
below.  The  intensity  I  is  given  directly  in  units  of  watts /cm^  in  our 
computer  program,  where  the  assumed  source  strength  is  1  watt/cm.  When 
I  is  multiplied  by  a  ,  the  absorption  coefficient,  the  result  oJ  is  the 
power  absorbed  per  unit  volume  ( watts/cm 3).  Since  I(a, a)  is  computed 
at  large  number  of  points  it  is  a  simple  matter  is  perform  the  surface 
integral  of  I  over  the  cross  section  of  the  rod,  which  yields  in 
watts/cm  the  power  absorbed  in  the  rod.  This  can  be  divided  by  the 
power  radiated,  N  watts/cm  from  all  pump  sources,  to  give  the  pumping 
efficiency .  The  formula  is 
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n  =  AfiJ*/  I  (a,  a)  ada  da 
_  2  I(a,a)  ada  da 


(41) 


This  method  of  calculating  n  is  not  as  simple  as  the  first  method,  but 
is  a  simple  addition  to  the  computer  program  that  calculates  intensity 
distributions .  For  comparison  we  calculated  the  n  for  the  last  case 
discussed  above  and  got  0.7335  using  Eq.  44  and  0.7492  using  Eq.  37. 
(So  far  we  have  not  found  the  cause  of  this  discrepancy.) 

Axial  Propagation. 

The  angle  6  that  the  ray  makes  with  the  transverse  plane  has 
so  far  been  assumed  to  be  zero.  Relaxing  this  restriction  makes  Figure 
2  a  projection  in  the  ( r, 8 )  plane  of  the  ray  in  three  dimensional 
(r,B,z)  space.  Thus  the  angles  8^.  and  8'.  are  defined  in  relation  to 
the  projections  of  the  rays  on  to  the  transverse  plane.  With  the  cylin¬ 
drical  symmetry  assumed  in  this  treatment,  a  ray  crossing  a  boundary 

between  media  and  will  make  angles  (90°  -  6^)  and  (90°  -  6 ^)  with 

the  cylinder  generator ;  and  phase  matching  along  this  generator  leads 
immediately  to 

n i  sin  & =  n^  sin  6,,  (42a) 

While  the  fact  that  the  two  rays  and  the  normal  are  coplanar  gives 

sin  8 j/sin  8^  =  tan  6^/tan  6^  (42b) 

The  ray  tracing  thus  gives: 

6^  =  constant  up  to  the  boundary,  and 

r  sin  8^  =  constant  up  to  the  boundary,  so  that  Bj  is  function  of  r. 
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Then  Eq.  42a  gives  a  new  & (called  5 )  and  Eq.  42b  gives  a  new 

Bfcalled  B  2)  •  Thence  6^  and  r  sinB^  remaan  constant  up  to  the  next 

boundary .  The  net  result  is  that  the  formula  for  each  8  .  and  B'.  in  Eq. 

JJ 

1  (as  used  in  Eq.  5 a  and  5b  for  example)  has  (n/nj)  replaced  by 
(tan  5  Jtan  &),  which  is  equivalent  to  the  substitution: 


—  -  (—)  cos  i/[l  -  (~)2  sin2  &  ]  *  =  M.(6)  ; 

n  .  n  .  n  .  7 

J  J  J 


where  6  is  defined  at  the  diode  array  source  along  with  $  as  the  angular 
coordinate  in  the  antenna  pattern.  (Recall  that  n  is  the  index  of 
refraction  of  the  outermost  medium.) 


In  propagating  from  a  diode  array  to  the  field  point  (a,u,z) 
inside  a  laser  rod,  the  ray  advances  axially  through  the  distance 


A z  =  z  -  z  .  =  I  (r  .  cos  8  .  -  r  .  ,  cos  B".  J  tan  &  . 

1  J  J  J'1  J’1  J 


where  the  last  r .  and  are  a  and  sin  (p/a)  respectively,  and 
where  the  last  term  will  also  have  a  "plus"  sign  rather  than  "minus"  in 
the  case  of  a  =  ( the  second  intersection  of  the  ray  with  the  cylin¬ 
der  r  =  a) .  Here  z ^  and  z  are  axial  coordinates  at  the  source  and  the 
field  point  respectively .  The  8 .  are  all  found  from  Eq.  1  with  sub¬ 
stitution  of  Eq.  43,  while  the  in  Eq.  44  are  given  by  Eq.  42a  as 


sin  6  .  =  (n/n  .)  sin  6  , 
J  J 


(42c) 


tan  5  .  =  M  .  tan  5 
J  J 


(42d) 
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It  is  necessary  in  the  three  dimensional  problem  to  find  both 
6  and  4>  that  identify  the  ray  through  the  field  point.  The  equations 
to  be  solved  for  (b,&)  are  first  Eqs .  5: 

a  =  b.  +  b  +  (B.  -  B'.)  -  y  ( 5a) 

1  i  j  j 


and 


where 


(*2  =  Oj  -  it  +  2y 


sin  S.  =  (nR  sin  b/n  .r  .)  M. 
J  J  J  J 


sin  B'.  =  (nR  sin  b/n.  ,  r  .)  M. 
J  J+1  J  J 


sin  y  =  (nR  sin  b/n^a)  M ^  =  p/a  ; 


(5b) 


and  then  Eq.  44: 

z  =  z.  +  l  (.  .  .)  M.  tan  6  .  (44) 

j  J 

Various  computational  methods  are  available  for  solving  for  (b, 6)  such 
as  Newton's  method,  generation  of  look-up  tables,  and  successive 
approximations  such  as  (b°,0),  (b°,f>°);  etc.  where 

(bn+*,6n)  solves  Eq.  5  with  6  fixed  at  &n  and  (bn,&n)  solves  Eq.  44 
with  4)  fixed  at  <l>n  , 

Taking  account  of  the  axial  component  of  propagation  will  have 
little  effect  on  the  estimates  of  pump  intensity  obtained  by  the  two- 
dimensional  propagation .  In  the  multipass  situation  involving  reflec¬ 
tion  of  the  transmitted  light  back  into  the  rod,  there  will  be  a  walkoff 
effect  that  will  mainly  smooth  out  the  intensity  variations  between 
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neighboring  arrays;  but  the  various  ray  angles  (and  associated  intensity 
calculations )  will  be  well  approximated  by  the  two-dimensional  theory. 
The  easiest  way  to  see  all  this  is  to  combine  Eqs .  42a  and  42b  to  give: 


sin  & 

sin  8  „ 
2 


n2  2  n2  2 

K~)  ♦  {(—) 

ni  ni 


(45a) 


which  for  small 


can  be  approximated  by 


sin  8^ 

7I^T2 


2 

tan 


(45b) 


This  will  be  less  than  a  one  percent  deviation  from  n2/n ^  for  the  range 
of  6  (a  few  degrees)  that  will  occur. 

Recall  that  the  two-dimensional  intensity  formula  (Eq.  10)  has 

on  its  right  side  terms  of  the  type  dW/dw  multiplied  by  the  transmission 

coefficient  T.  The  numerator  may  be  considered  to  be  (dW/d<b)d$  , 

representing  the  power  radiated  by  the  source  i  into  the  angle  d< J>  ; 

while  the  denominator  similarly  is  (a  -  p  )  |  c?ac/cf<J»  |  c?4>  ,  representing 

the  line  element  perpendicular  to  the  ray  through  which  the  power  dW 

flows.  In  three  dimensions  the  formula  is  readily  patched  by  changing 

2 

the  numerator  to  (3  W/3$d6) d<t>d6  and  putting  the  additional  factor 
(dz/dS)  cos  6jd6  in  the  demoninator  in  order  to  generate  the  area 
element  through  which  power  dW  flows.  The  net  result  is: 


T i (dWj/dfy) (dW2/d&) 

I(a,a, z)  =  L  - 

1  {a2  -  (M sin  4>)2}  3  ||^|  {I  -  (n  sin  S/n ^)2}  3 

(46a) 
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Here  we  have  replaced  the  second  derivative  of  W  in  the  numerator  by  the 
product  of  the  two  first  derivatives  on  the  assumption  that  &($,&) 
factorizes ,  which  is  almost  certainly  the  case  with  diode  arrays .  Since 
there  is  not  much  depolarization  for  small  &  the  exact  ($,(>)  -  depend¬ 
ence  of  T ^  are  not  shown.  The  values  given  by  Eq.  11  and  12  are  close 
enough . 


Eq.  46a  is  good  for  a  source  having  zero  extent  in  the  z- 
direction,  and  needs  to  be  modified  for  a  diode  array.  The  z-coordinate 
in  the  array  is  written  z .  as  in  Eq.  44.  If  the  array  is  coherent  the 
radiation  at  the  field  point  will  be  of  near- field  character,  and  the 
factor  dW^/d6  in  the  numerator  combines  with  f3/36 ) (z-z  .)  in  the  denom- 

2  l 

inator  to  give  dW^f&f/dz^  ,  which  characterizes  the  source  radiation  at 
direction  6  .  (In  Equation  46a,  z ^  was  assumed  to  be  zero.)  A  given 
field  point  may  be  illuminated  by  two  separate  coherent  arrays  if  the 
arrays  each  radiate  at  two  distinct  angles  (e.g.  ±  6 )  . 

If  on  the  other  hand  the  array  is  incoherent  one  must  add  the 
illuminations  from  different  elements  along  the  array,  and  form  a  sum 
over  z ^  of  terms  of  the  type  that  contain  5  in  Eq.  46a.  Expressing  the 
result  as  an  integral  gives  the  factor 


OW^/dz^  *  (dW2/d6)  |  36/3zi  |  } 
{ 1  -  (n  sin  5 /n^}  1 


dz 


(46b) 


where  6  is  given  as  a  function  of  z for  fixed  The  factors  contain¬ 

ing  4>  in  Eq .  46a  remain  the  same. 
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Convergent  Pump  Beams  and  Finite  Apertures. 


The  analysis  referred  to  in  Fig.  1  needs  to  be  extended  to 
cover  the  case  where  the  light  from  each  diode  array  is  made  weakly  con¬ 
vergent  by  passing  through  a  positive  lens.  The  resulting  beam  passes 

through  an  aperti  re  in  the  reflecting  surface  at  radius  R.  (Fig.  1  is 

still  appropriate  for  a  diverging  beam;  but  if  the  reflecting  surface  is 
located  at  a  radius  less  than  that  locating  the  virtual  sources,  the 

reflector  will  have  holes  in  it,  and  the  present  analysis  can  be  used  to 
treat  this  case  as  well.)  Fig.  14  shows  the  aperture  A^AA2  in  the 

reflector ,  with  pump  light  converging  towards  a  point  C  on  the  diameter 
at  the  distance  D  from  the  entry  point  A.  Angular  positions  on  the 
wavefront  are  labeled  by  \  referred  to  C  (towards  which  the  rays  are 
directed) ,  while  the  aperture  itself  is  labeled  by  the  variable 
4k,  with  4k  equal  to  4kq  at  the  center.  The  index  i  counts  the  pump 

sources.  4> .  is  the  same  quantity  as  0,,  in  Eq.  5a.  As  before  4> (i) 

1  n 

abbreviated  to  4>  ,  denotes  the  angle  between  the  ray  and  the  radius,  and 

in  accordance  with  the  sign  convention  chosen,  is  positive  as  shown  in 

the  figure.  There  is  a  one-to-one  correspondence  between  4>  and 

f4>.-4>.  ),  and  between  4>  and  y  •  The  offset  of  the  ray  projection  from 
1  10 

the  center  0  is 


R  sin  4>  =  (D  -  R)  sin  x 


while 


4> .  -  <b  .  =  y  + 

vio  * 


.  -I  ,R  sin  4>, 
♦  +  S1*  l~D~~R} 


(47) 

(48) 

(49) 
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APERTURE 


WAVEFRONT 


REFLECTOR 


OA  =  R 
AC  =  D 
OC  =  R  -  D 
R  sin  0  *  (R  -  D)  sin  x 

0i  -  0io  =  X  +  0 


C88-8955 

Figure  14.  Geometry  of  wavefront  converging  on  C  passing  through 
aperture  AjAAg  in  reflector  centered  on  0.  $  labels 
the  point  in  the  aperture  at  which  the  ray  direction  is  <t> 
(shown  positive),  x  labels  the  point  on  the  wavefront  in 
terms  of  which  the  antenna  pattern  of  the  source  is 
defined. 


50 


These  relations  hold  irrespective  of  the  sign  of  D  or  its  size  in 
relation  to  R.  The  interesting  case  for  D  positive  is  with  D  large 

compared  to  R  so  that  no  focus  forms  in  the  rod.  Thus  from  Eq.  49, 
d<t>^/d$  is  positive. 

Notice  also  that  the  pump  source  distribution  expressed  in 
terms  of  x  is  readily  transformed  to  the  variable  $  : 

dW  =  (dW/dx)  d\  =  (dw/dx)  (d\/d<b)d$  -  (dW/d$) d$ 

so  that  Eq.  10  for  the  intensity  at  a  point  in  the  rod  is  still  valid. 
For  the  small  values  of  $  and  x  involved,  Eqs .  47  and  49  are  essentially 
linear  relations  amongst  the  variables . 

The  method  for  finding  the  value  of  corresponding  to  a  for  a 
given  radius  a  is  now  a  little  different .  Previously,  given  a  ,  we 
imagine  a  horizontal  line  at  height  a  on  Figure  4  intersecting  the 
curve  a(<f>)  to  identify  $  .  The  algebraic  problem  to  solve  was  Eq.  7: 

a($)  -a  -  0 

The  curves  af<t>J  in  Fig.  4  are  drawn  for  $ ^  =  0  ,  and  have  *<t >  as  the 
horizontal  axis.  They  are  identical  to  plots  of  ($ .  -  a)  versus  -4>  . 
Now  the  problem  is  to  find  that  solves  Eqs.  7  and  49  simultaneously, 
i .  e. 

4k  -  a  =  <jK  -  a($)  (7a) 

and  4> .  —  4> .  =  $  +  sin  *  (44) 

i  10  D  -  R 

The  method  of  solution  is  illustrated  graphically  in  Fig.  15, 
from  which  it  is  clear  that  the  problem  to  be  solved  is  to  find  (fit hat 
satisfies  the  following  equation: 
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Figure  15.  Graphical  illustration  of  the  simultaneous  solution  of 
Equations  7a  and  49.  a,  b,  c  illustrate  Eqs.  5,  7a,  and 
44  respectively.  The  solution  .)  is  given  by  the 
intersection  in  d.  1 
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(50) 


0  =  [<J>  -  a($)  ]  *  a  -  (<fr  *  4)  +  sin  1  {—  ' Sin|}  ] 

l  lO  D  -  R 

The  quantity  in  the  first  square  brackets  is  the  function  obtained  from 

Eq .  5  by  moving  $  .  (which  is  0  )  to  the  left  side  of  the  equation  and 
l  M 

then  reversing  the  sign.  The  next  term  a  is  a  number  giving  the  angular 
coordinate  of  the  field  point.  The  quantity  in  the  second  square 
brackets  is  the  function  identified  by  Eq.  49.  In  shorter 

notation,  Eq.  50  is 

0  =  Z  &'.(*)  -  Z  B  .(*)  *  y(t)  +  a  -  *  .(*)  (50a) 

,•  J  •  J  -i 


Here  the  S',  and  8.  are  identified  in  Eqs.  1  and  5, 
J  J 

to  last  term  of  Eq.  5a. 


and  -y  is  the  second 


After  solving  Eq.  50  for  $  it  is  an  easy  matter  to  apply  Eq. 
10  for  the  intensity .  The  terms  in  the  numerator  are  now:  T(dW/d\) 
(d\/d$)  (where  the  subscripts  i  are  dropped  for  simplicity).  The 
denominator  is  now  formally  the  same,  but  in  the  computation  it  is 
necessary  to  add  to  the  formerly  used  expression  for  da/d$  the  new  term: 


d<t»i/d<|)  =  1  *  dx/d$  (51) 

where  we  have  used  Eq.  48.  A  simple  check  on  the  theory  is  to  treat  the 
case  shown  in  Figure  16,  even  though  the  pump  waves  are  diverging  in 
this  case,  because  we  already  have  the  results  in  terms  of  the  earlier 
analysis  wherein  C  is  located  on  the  circle  of  radius  R'  .  Let  -f(to) 
represent  the  functional  form  in  the  first  square  brackets  in  Eq.  44; 
i.e.  f  is  the  expression  given  by  Eq.  5.  Then  the  modified  form  of  Eq. 
10  contains 
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Figure  16.  Geometry  of  aperture  with  diverging  wave,  where  D  is 

negative.  The  radius  of  the  reflector  is  less  than  the 
radius  of  the  surface  containing  the  sources.  Here  r' 
replaces  R  in  Eq.  1  and  other  equations  in  that  section. 
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( dW/dy)  (dx/d<j>)  _  dW/d x _ 

(d/d<$>)(f  +  (d$/dx)  (d/d$)  (f  +  4^,1 

dW/dx 

da/dx 

thus  yielding  the  correct  result  referred  to  angles  x  at  the  source,  as 
expected .  The  significance  of  this  result  is  trivial  and  obvious  when  D 
is  a  small  positive  quantity  (smaller  than  R  -  r^  in  the  case  of  four 
concentric  media),  or  any  arbitrary  negative  quantity,  because  it  means 
that  we  simply  deal  with  the  positions  and  antenna  patterns  of  the  pump 
diode  images  in  all  analysis .  The  real  significance  is  that  for  large  D 
(convergent  beam  from  the  source-lens  combinations)  the  original  for¬ 
mulae  are  all  still  good  although  R  is  assigned  a  negative  value  in  all 
computations  i.e.  R  is  replaced  by  (R  -  D) ,  which  is  now  negative,  and 
this  new  parameter  is  renamed  R;  and  x  is  renamed  0  .  Thus  in  most 
computations,  where  lens  aberrations  are  unimportant,  it  will  be 
unnecessary  to  deal  with  Eqs.  49-51. 
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SBCTION  2 


CAUSTICS 


The  results  of  Section  1  can  be  applied  readily  to  calculate 
the  pumping  intensity  disribution  inside  a  cylindrical  rod  very  quickly, 
and  probably  quite  accurately  except  in  the  neighborhood  of  the 
caustics .  This  section  presents  the  elementary  geometrical  theory  for 
predicting  the  location  of  the  caustics,  develops  the  wave  theory  for 
calculating  the  intensities  in  the  regions  where  the  geometrical  theory 
fails,  and  discusses  measures  that  reduce  the  anomalous  intensities 
associated  with  the  caustics . 

Location  of  Caustics 


A  caustic  is  the  locus  of  the  point  of  intersection  of 
neighboring  rays,  which  are  tangential  to  the  caustic  curve.  In  our 
notation  the  caustic  is  characterized  by  da^/d^  being  equal  to  zero  at 
fixed  radius  a.  The  existence  of  the  caustics  is  confirmed  immediately 
by  plotting  a at  the  outer  surface  of  the  rod  i.e.  at  a  =  r^,  as 
illustrated  in  Figure  17.  The  minimum  of  the  curve  shown  locates  the 
point  where  the  caustic  meets  the  surface.  Minima  on  similar  curves  for 
successively  smaller  values  of  a  determine  polar  coordinates  along  the 
caustic.  Figure  18  shows  the  resutant  caustic  curves.  The  values  of 
the  coordinates  for  the  upper  curve  are  tabulated  together  with  corres¬ 
ponding  values  of  6  in  Table  III.  The  caustics  can  occur  only  for 

points  (a, a)  where  a  exceeds  p  ,  because  there  the  a($)  curve  has  two 

M 

branches  as  in  Figure  4.  The  single  branch  curves  of  Figure  4  have  no 
maximum  or  minimum  value  of  a  other  than  0  and  180  degrees. 
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Caustic  cu 
parameters 
efficiency 
also  drawn 


TABLE  III 


Coordinates  along  the  upper  caustic  of  Figure  18 


a 

(X 

<t> 

( cm) 

( deg. ) 

(deg. 

0. 3175 

164.03 

5.0 

0.  30 

161.43 

5.1 

0.29 

159.69 

5.2 

0.28 

157.69 

5.3 

0.27 

155.36 

5.4 

0.26 

152.59 

5.525 

0.25 

149.13 

5.645 

0.24 

144.40 

5.80  5 

0.235 

140. 85 

5.93 
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Caustics  are  to  be  expected  in  cases  of  practical  interest  in 
this  work,  and  will  not  be  eliminated  by  the  addition  of  other  concen¬ 
tric  media  to  the  simple  geometry  depicted  in  Figure  1.  Referring  to 
Eg.  8b  for  do  /d$  it  is  seen  that  as  $  approaches  $  (which  is 
sin'l (n2r ^/nR) )  the  second  positive  term  within  the  braces  approaches 
infinity  as  illustrated  in  Figs.  4  and  17,  for  which  a  =  p  .  There- 
fore,  unless  da^/db  is  positive  at  $  equal  to  zero,  it  will  be  zero 

somewhere  in  the  interval  0  <  $  <  $  .  Substituting  <J>  =  0  in  Eq.  8  with 

Ml 

a  equal  to  r ^  shows  that  the  condition  for  no  caustic  is  that 


2 


nlrl 


for  four  media,  with  obvious  extension  to  larger  numbers  of  media. 

Replacing  vacuum  in  any  region  r'  <  r  <  r'  by  medium  with  index  n'  adds 

the  negative  quantity  -  1)(~  -  — 7)  to  the  right  hand  side  of  the 

n  rj,  r  2 

inequality ;  and  therefore  the  maxxmim  value  that  the  right  hand  side  can 

have  is  f—  -  j-)  ,  with  limiting  value  1/r Comparison  with  the  left 
rl  R 

side  shows  that  we  would  need  to  have  n^  >  2  ,  which  occurs  for  the 
niobates  amongst  the  common  laser  hosts,  and  maybe  something  else.  For 
YAG  or  YLF  or  anything  else  with  n ^  <  2,  caustics  will  occur  with  the 
concentric  geometry .  The  only  ways  to  keep  the  caustics  from  forming  in 
the  gain  region  are  to  limit  the  range  of  angles  from  the  pump  source  by 
suitable  design  of  the  coupling  optics,  or  to  have  the  laser  rod 
material  at  radii  exceeding  (n^r^/n^)  undoped. 


Because  the  upper  branch  of  the  a  versus  <t>  curve  has  the 

infinite  slope  at  its  extremity  the  inner  end  of  the  caustic  is 

identified  with  the  point  at  the  angular  break  on  curve  4  of  Fig.  4, 

which  is  at  the  mid-point  of  the  chord  marking  the  critical  ray  (i.e. 

$  =  <j>.J  •  Each  pump  source  generates  two  caustics  symmetrically  about 
n 

its  central  ray.  Caustics  from  different  pump  sources  intersect, 
doubling  locally  the  intensity  anomaly.  (With  five  pump  sources  the 
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angular  positions  of  the  intersections  are  zero  and  multiples  of  72 
deg.)  We  have  no  mathematical  formula  that  generates  the  caustic 
curves,  but  it  is  a  simple  matter  to  construct  a  table  of  coordinates 
like  Table  III  once  a  calculator  has  been  programmed  to  generate 
ot^(<f;a)  .  The  table  then  serves  as  a  guide  for  selecting  a  cut  through 
the  geometrical  optics  caustic  along  which  intensity  can  be  computed . 

In  terms  of  the  multipass  analysis  used  in  Section  1  to 
account  for  reflection  of  the  transmitted  light  back  into  the  laser  rod, 
it  should  be  noted  that  caustics  are  formed  by  the  initial  set  of  rays 
(before  reflection)  only.  The  multipass  formula  for  intensity,  which  is 
believed  to  be  a  good  approximation  in  general ,  is  not  used  for  investi¬ 
gation  of  intensity  distribution  near  caustics . 

WAVE  THEORY  OP  THE  CAUSTIC 

The  Diffraction  Integral. 

The  geometrical  optics  (G.O.)  theory  identified  the  caustic  as 
the  locus  of  points  (a, a)  where  3a/3$  is  zero,  because  this  derivative 
appears  in  the  denominator  of  the  expression  for  the  intensity .  The 
wave  theory  works  by  expressing  the  field  amplitude  at  a  point  in  the 
region  of  the  caustic  in  terms  of  a  diffraction  integral.  The  surface 
of  integration  is  any  surface  over  which  the  field  is  defined  accurately 
by  the  G.O.  theory.  The  surface  chosen  is  the  inside  surface  of  the 
laser  rod  at  radius  r ^ .  The  complex  amplitudes  at  points  (r  ,Q,z)  on 
the  surface  are  given  by 

iks ' 

(ri,e,z)  =  [I(r ltQ,z)]  ‘  e  2  (101) 


wherein  the  units  of  A  are  such  that  the  square  of  the  magnitude  gives 
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the  intensity  inside  the  medium  of  refractive  index  n^.  The  optical 
path  S'  from  the  source  to  the  rod  surface  is  given  by 

S' 2  =  S'('ri,  Q,z)  =  J  nds  =  I  n^.  s (102) 

evaluated  along  the  ray  AQ.  The  geometry  is  shown  in  Figure  19  where  A 
denotes  the  point  or  line  source,  Q  the  point  (r^,0,z)  that  generates 
the  Huygens  wavelet,  and  P  the  field  point  (a,a,0) 

The  field  at  P  is  then 


(ikn  -a/2) S' 

A(P)  =  (-in2/\)  J  J  A(Q)  e  1  0(B,z)dq/S'  (103) 

where  dq  is  the  area  element  at  Q,  S'  is  the  distance  QP,  and  0(0)  is 
the  obliquity  factor  for  the  Huygens  wave  source  at  Q: 

0(0, z)  =  j  (1  +  cos  t)cos  (104) 

where  z(0,z)  is  the  angle  between  QP  and  the  refracted  G.O.  ray  at  Q, 
and  8^  is  the  angle  of  refraction  inside  the  medium  n^,  i.e.  the  angle 
between  the  surface  area  element  and  the  wavefront  element. 

Fig.  19  is  projection  on  the  transverse  plane  and  thus  does 
not  show  the  vertical  dimension  z.  If  S'  is  the  transverse  component  of 
S'  .  then 


(S') 


2 


2 

z 


and  for  the  small  range  of  z  involved  in  the  integration,  the  approxi¬ 
mation 
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Figure  19.  Projection  of  ray  path  on  plane  perpendicular  to  z-axis. 

The  shaded  area  is  the  surface  of  integration  used  in  the 
diffraction  calculation.  P  is  the  field  point,  and  A  the 
source.  The  path  from  A  to  Q  is  determined  by  geometrical 
optics.  QP  is  the  line  along  which  the  Huygens  wavelets 
propagate. 


(105a) 


S'  =  Si  +  z  / 2S z 


is  made.  Similarly 


S'  =  S2  +  z  / 2L 


(105b) 


where  L  has  the  dimension  of  length.  The  quantity  L  is  infinite  if  A  is 
a  line  source ,  i.e.  A(Q )  then  has  as  its  phase  factor  exp( ikS ) 
Finally,  with  the  area  element  expressed  as 


dq  -  r^dQ  dz 


the  diffraction  integral  becomes: 


ik(S  *  n  S  )  -cS/2 

l (P)  =  (-in  /\)r  f  dQ  •  e  * 


I  44  7  7  1 

J  dz  •  0(d,z)(S2  +  zZ  /  2S  x )  1  [Ij  (ty  8,  z)  ]2  exp[  {i  —  [  —  +  ~)  -  ~"}z  1 


(106) 


The  phase  in  the  second  integral  is  stationary  at  z  -  0  and  varies 
rapidly  on  either  side  of  this  point,  so  that  this  integral  is  approxi¬ 
mated  closely  by  the  expression 


J  .  2  2 

Integal  =  [I  fr  ,9,0)]  T  0(d,0)/S7)  f  "  elK  z  dz 
11  l  -co 


(107a) 


and  this  last  integral  is  equal  to  (in)  1  /K  where 
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K  =  K(Q)  =  t (k/2)(n1/S1  +  1/L)  +  ia/4S1\ 


(107b) 


(The  term  proportional  to  a  is  neglible  compared  to  the  term  in  k.) 

A 

Because  of  this  last  expression  one  gets  the  S ^  factor  in  the  denomi¬ 
nator  of  the  0  integral  in  Eg.  106  for  cylindrical  waves,  i.e.  when  1/L 
is  zero.  At  this  point  we  specialize  the  calculation  to  this  case  and 
write  for  Eq.  107: 

Integral  =  [I1(rJL,Q,0)]  *  0(6,0)  (107c) 

Substitution  of  Eg.  107  into  Eq.  106  gives 


3vi/4 
e  r  ,n 


A(P)  = 


iii(Q ) 

11  J  A(Q)  e  ^|oye  0(9.0) 


d  9 


(108a) 


( n 2\) 


A 

S1 


for  the  field  amplitude,  and 


I(P) 


2 

rini 


\J(P)\2 


(108b) 


for  the  intensity,  where  J(P)  is  the  integral  in  Eg.  108a.  Note  that  in 

A 

the  expression  for  J(P)  the  amplitude  [1^  (r^  9, 0)  ]  *  has  been  written 
A(Q).  The  phase  i>(9)  is  given  by 


i>(8)/k  =  S  ^  +  nS 


(109) 


TTC-1529-R 


65 


for  large  enough  k  the  phase  will  change  so  rapidly  as  a  function  of  0 
(away  from  its  stationary  point )  that  again  the  stationary  phase 
approximation  can  be  used  and  the  whole  behaviour  of  I(P)  is  determined 
by  the  way  that  'V(Q)  varies  in  form  from  one  point  P  to  another.  In  the 
case  of  practical  interest  here,  k  is  not  large  enough  for  this  approxi¬ 
mation,  but  it  is  fruitful  to  investigate  this  asympotic  form  because  it 
can  give  a  quick  approximation  to  the  integral. 

The  ratio  of  intensities  at  two  field  points  is  the  ratio  of 
the  corresponding  quantities  given  by  Eq .  108b  and  it  is  convenient  to 
compare  intensities  to  those  on  the  axis  of  the  rod  where  the  C.O. 
theory  gives  an  accurate  estimate . 

The  Phase  Factor. 

The  variation  of  't/k  with  0  is  shown  graphically  in  Figure  20 
for  four  field  points,  two  in  the  illuminated  zone,  one  on  the  G.O. 
caustic,  and  one  in  the  G.O.  shadow  respectively,  reading  from  top  to 
bottom.  (All  the  curves  have  a  cut-off  at  the  right  hand  end  because  of 
the  critical  ray.)  The  top  curve  has  a  minimum  that  identifies  0  for 
the  G.O.  ray.  The  second  curve  has  two  points  of  stationary  phase  cor¬ 
responding  to  the  maximum  and  minimum,  thus  identifying  two  rays  through 
the  field  point.  One  ray  satisfies  Fermat's  principle  with  a  minimum 
and  the  other  satisfies  it  with  a  maximum.  Interference  occurs  because 
of  the  two  rays.  Each  of  the  curves  can  be  represented  by  an  equation 
that  contains  primarily  linear  and  cubic  terms  when  expanded  about  a 
8-value  of  a  little  less  than  one  radian.  On  the  third  curve,  for  a 
point  located  on  the  G.O.  caustic,  the  linear  term  has  zero  coeffici¬ 
ent.  The  bottom  curve,  which  has  linear  and  cubic  terms  of  the  same 
sign,  has  no  extremum  and  therefore  no  G.O.  ray. 
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4.9 


a  =  0.27  CM 


88-8960 


Figure  20.  Plot  of  optical  path  length  defined  by  Eq.  109  versus  e  , 
where  is  shown  in  Fig.  19  and  S2  is  the  optical  path 
AQ.  The  geometry  used  in  this  illustration  has  four 
regions  with  radii  0.3275,  0.5,  0.65  and  4  cm  and  indexes 
of  refraction  1.823,  1.333,  1.5  and  1.  The  solid  triangles 
point  to  the  extrema  on  the  curves. 


In  the  neighborhood  of  the  caustic  the  field  intensity  will  be 
greatest  at  the  point  for  which  the  phase  curve  has  the  largest  range  of 
0  that  limits  the  range  in  <li  to  approximately  it  .  The  phase  will  thus 
have  the  form  of  the  second  curve  of  Figure  20  with  a  vertical  offset 
between  maximum  and  minimum  of  somewhat  less  than  it  ,  and  thus  the 
corresponding  point  lies  not  exactly  on  the  G.O.  caustic  but  just  inside 
the  illuminated  zone.  This  is  analogous  to  finding  that  the  maximum  of 
the  quantum  mechanical  wave  function  for  a  particle  being  turned  on  a 
potential  slope  occurs  not  at  the  classical  turning  point  ( where  the 
wavef unction  has  an  inflexion  point )  but  on  the  downhill  side  of  this 
point . 


If  the  amplitude  of  the  Huygens  wavelets  were  independent  of 
0  and  the  other  factors  in  the  integrand  of  Eq.  108a  constant,  the 
integrand  would  be  of  the  form  exp  i(bQ  +  a 3  9^ )  which  yields  Airy 
functions  ( Bessel  functions  of  order  one-third)  when  integrated  over  a 
large  number  of  cycles.  The  ratio  h/a  varies  from  one  curve  to  the  next 
in  Fig.  20.  In  the  shadow  region  were  b  has  the  same  sign  as  a  one 

would  obtain  the  Bessel  function  )>  while  in  the  illuminated  zone 

where  b/a  is  negative  the  combination  +  J  1/3^^  would  occur. 

In  each  case  the  argument  y  is  (\2b/3a\ (\b/3a\)  1  }  .  The  integral  is 

readily  obtained  by  combining  Formulae  1  and  2  of  I.S.  Gradshteyn  and 
I.H.  Ryzhik  "Table  of  Integrals,  series,  and  Products"  section  3.695. 
(Academic  Press  1980). 


Hence  in  the  limit  of  large  k  the  intensity  at  P  as  given  by 
Eqs .  108  and  the  referenced  formulae  is 

-aS 

I(P)  =  fr^/X  S 1 )  •  | 4tt 2b/27a3\  •  {J  fy )  +  J-l/3(y)}?  '  I(Q)  e  1  X 
2 

<  cos  S 1(Q)> 

(110) 
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in  the  illuminated  zone,  with  a  similar  formula  having  (y)}  for 

the  shadow  region.  Here  <  .  .  >  denotes  an  average  value  of  the  cosine- 
squared  factor  e.g.  its  value  at  a  point  midway  between  the  maximum  and 
minimum  phase  points. 


Estimate  of  the  Peak  Intensity. 

Although  modern  optics  texts  do  not  treat  the  subject  of 
caustics  it  is  no  accident  that  the  Airy  function  is  tailored  for  their 
description.  The  Airy  function  is  defined  by 


Ai(-x)  =  (i/3)  x  1  { J1/3(y )  + 


(Ilia) 


Ai(x)  =  tt  (x/3)  1  Kj/3(y^  ' 


(111b) 


where  y  =  (2/3)  x2^2.  (See  H.  Abramowitz  and  I. A.  Stegun  (Eds.)  "Hand¬ 
book  of  Mathematical  Functions,  "  pp.  446,  447,  478;  especially  Table 

2  2 

10.13.).  It  is  not  related  to  the  Airy  pattern  [J^( x)]  /x  describing 
diffraction  by  a  circular  aperture. 


The  expression  111a  will  be  used  here,  wherein  y  is 
3 /2 

2(-b/3a)  ,  which  is  positive  because  b  is  positive  and  a  negative  in 

the  region  of  interest.  In  terms  of  Eg.  111a, 


c  ao  3  3  1/3  1/3 

J  exp  i(a  x  +  bx)dx  =  (~2v/3  a)  Ai(b/3  a)  (112) 

*  — a> 

so  that  Eq.  110  now  reads: 


I(P)  =  (r^n^XS^  (4*2/32/3a2)[Ai(b/31/3a)  ]2  I(Q)  e  1  cos2B1(0 )  (113) 
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The  maximum  value  of  Ai(-x)  occurs  at  its  first  maximum,  which 
is  0.53565666  at  x  =  1.01879297.  Putting  this  value  into  Eq.  13  gives 

2  2  2 

HP)  (Maximum)  =  (5.4457/a  ifr^n^/XS^  I(Q)  e  cos  &1(Q)  (114) 

The  specific  curve  of  this  type  shown  in  Figure  20  that  corresponds  to 
(-b/31^3 a)  being  equal  to  x  above  has  a  vertical  displacement  between 
maximum  and  minimum  values  (equal  to  5(b/3a)^^ )  of  1.7139  radians 
(98.2  degrees),  thus  putting  the  maximum  just  inside  the  illuminated 
region  as  expected . 

This  estimate  of  the  maximum  intensity  expected  in  the  caustic 
region  has  as  its  principal  assumption  the  constancy  of  the  intensity 
I (Q)  on  the  inside  surface  of  the  rod  over  a  region  that  takes  in  a 
large  multiple  of  2 ir  in  the  phase  function  ty(Q)  on  each  side  of  its 
slowly  varying  part,  i.e.  that  there  be  a  large  number  of  Fresnel  zones 
on  each  side  of  the  pair  of  G.O.  rays.  The  assumption  is  not  neces¬ 
sarily  a  good  one  in  cases  of  practical  interest.  The  various  curves  in 
Fig.  20  vary  principally  in  their  values  of  b.  The  coefficient  a  varies 
more  slowly  and  the  value  to  be  used  in  Eq.  114  can  be  taken  from  the 
third  curve  i.e.  the  one  with  b  equal  to  zero.  The  coefficient  a 3  is 
determined  from  the  table  of  values  used  to  plot  ty/k  . 


Fig.  20  refers  to  points  close  to  the  caustic  of  Fig.  18  in 
the  previous  section  "Transverse  Pumping  -  Geometrical  Theory. "  There, 

X  =  0.8  micron,  a 3  =  3108,  r  =  .3175,  n  =  1.823  and  we  get  from  Eq. 
114 


I(P) (Maximum)  =  141.5  I(Q) 


(Q) 


For  the  other  parameters  applying  to  this  particular  configuration, 
including  a  =  3  cm  * ,  and  6^  =  40.57  deg.  ,  the  above  relation  becomes 
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I(P)  (Max. )  =39.2  1(0) 


in  terms  of  the  intensity  on  the  axis,  for  the  caustic  at  radius  0.27 
cm.  As  will  be  seen  below  this  last  expression  over  estimates  the 
intensity  by  about  a  factor  of  4. 

Computation  of  the  Integral . 


A  computer  program  vas  written  by  our  subcontractor  (Pharos 
Technical  Enterprises)  to  calculate  the  intensity  in  caustics  of  the 
pump  field  on  IBM  PC  computers  and  clones.  The  program  is  written  in 
Microsoft  QuickBasic  4.00b  and  is  called  8PGM2100 .  727 .  A  stand-alone 
executable  version  is  available  which  runs  outside  the  QuickBasic 
environment,  called  PG88-007.EXE.  The  stand-alone  version  runs  approxi¬ 
mately  twice  as  fast  as  the  QuickBasic  version.  Either  runs  with  or 
without  an  8087  (or  other)  numeric  coprocessor  with  no  change  in 
accuracy.  Software  emulation  of  the  8087  's  IEEE  floating-point 
arithmetic  format  is  done  in  the  absence  of  the  hardware  numeric 
coprocessor,  slowing  the  program  by  a  factor  of  approximately  30. 

The  program  calculates  the  Fresnel  diffraction  integral  for 
two  dimensional  pumping  geometries .  The  program  does  not  (but  could 
easily  be  modified  to)  dynamically  adjust  its  step  size.  As  a  result  it 
may  be  used  only  for  problems  closely  approximating  the  case  for  which 
it  was  written.  The  program  calculates  the  pump  intensity  at  the  field 
point  (r,8)  within  the  laser  rod  for  a  single  diode  located  at  a  radius 
approximately  6  times  the  radius  of  the  cooling  jacket  away.  The  number 
of  different  media  in  the  pump  path  is  arbitrary ;  4  were  used  in  the 
test  cases  (air,  water  jacket,  coolant,  laser  rod).  The  intensity  is 
calculated  for  the  first  pass  of  the  pump  light  through  the  rod. 
Similarly,  it  is  assumed  that  pump  fields  from  multiple  diodes  can  be 
added  using  geometrical  optics  approximations . 
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The  program  divides  the  problem  into  4  parts:  data  input, 
geometrical  optics  calculations,  diffraction  calculations ,  and  data  out¬ 
put  . 


Data  input  may  be  from  a  file  or  from  the  keyboard .  The 
system  may  be  edited  and/or  saved  in  a  file  for  reuse.  The  program 
handles  routine  errors  in  file  usage,  such  as  specification  of  non- 
existant  directory  or  file  names. 

The  method  of  the  program  is  to  select  a  radius  from  the 
center  of  the  rod  through  the  diode.  The  far-field  emission  pattern  of 
the  diode  is  then  projected  from  the  diode,  through  the  intervening 
media  ( if  any),  to  the  inner  surface  of  the  rod.  The  pump  field  is 
calculated  at  discrete,  fixed-step  intervals  around  the  periphery  of  the 
rod.  (This  needs  to  be  generalized  to  dynamically  optimize  the  step 
size.)  This  step  is  called  the  geometrical  factors  calculation,  and 
requires  about  1  to  2  minutes  for  a  4  media  geometry  on  a  PC.  It  uses 
the  G.O.  theory  of  Section  1.  The  E  field  is  presumed  to  be  parallel  to 
the  axis  of  the  system,  which  corresponds  to  the  normal  polarization  of 
diode  arrays  oriented  so  that  the  long  axis  of  the  chips  are  parallel 
with  the  system  axis. 

The  user  may  specify  a  range  of  field  points,  in  either  or 
both  angle  and  radius,  for  the  program  to  sept  through  automatically . 
For  each  field  point  the  program  integrates  the  contributions  from  all 
illuminated  portions  of  the  rod  periphery ,  using  linear  interpolation 
between  the  points  along  the  periphery .  The  use  of  linear  interpolation 
increases  the  number  of  points  along  the  periphery  that  must  be  cal¬ 
culated  in  the  geometrical  portion  of  the  program,  but  minimizes  the 
computation  time  in  the  diffraction  portion  of  the  program,  striving  for 
a  good  balance.  The  integration  is  done  in  fixed  steps  of  angle  around 
the  rod;  this  needs  to  be  generalized  so  that  the  program  can  optimize 
its  computation  time  and  handle  a  large  range  of  cases.  The  diffraction 
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integral  includes  all  the  necessary  terms  to  account  for  the  cylindrical 
geometry,  and  includes  the  obliquity  factor,  which  is  a  small  (of  order 
10~6)  adjustment  in  the  cases  tested.  It  does  not  include  a  term  to 
account  for  the  fact  that  H  vectors  are  not  all  parallel  at  the  field 
point . 


The  data  output  is  to  the  screen,  to  a  printer  (if  specified) , 
and  to  two  disk  files  (if  specified) .  Two  output  files  are  used  to  pro¬ 
tect  from  data  loss  during  unattended  operation  should  a  power  failure 
occur  during  the  short  period  that  data  are  being  written  to  a  file. 

In  order  to  display  the  data,  a  group  of  complementary  pro¬ 
grams  ( 8PGM0845 .  716 ,  et  al)  were  written  in  HathCad  2.0 6.  They  directly 
read  the  data  files  written  by  the  diffraction  integral  program. 
HathCad  provides  a  variety  of  tools  for  analyzing  and  plotting  the  data, 
and  there  seemed  no  need  to  provide  an  inadequate  subset  of  its  capa¬ 
bilities  in  the  diffraction  integral  program.  The  files  can  also  be 
read  into  other  analysis  and  plotting  programs,  spread  sheets,  or  word 
processor  programs ,  if  desired. 

The  diffraction  integral  program  has  been  extensively  cross 
checked  against  geometrical  optics  calculations ,  both  outside  of  the  rod 
and  within  the  rod.  The  agreement  has  been  found  to  be  more  than  ade¬ 
quate,  of  the  order  of  10~4  or  better  for  the  cases  tested.  The  penalty 
for  this  accuracy ,  and  the  lack  of  dynamic  step-size  adjustment  is 
slower  computation .  The  diffraction  and  output  portions  of  the  program 
require  3.5  to  7  minutes  per  field  point  on  a  PC/XT/AT  class  machine 
with  a  numeric  coprocessor .  The  geometrical  optics  programs  used  for 
comparison  were  written  by  another  programer,  in  a  different  version  of 
BASIC  (interpreted  BASICA) ,  utilize  a  different  arithmetic  (Microsoft 
Binary  Format  instead  of  IEEE  format),  and  run  on  a  different  machine. 
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Elsewhere  in  this  report  the  geometrical  propagation  equations 
and  the  form  of  the  Fresnel  diffraction  integral  used  are  discussed. 
There  are  also  displays  of  typical  data  including  comparisons  to  the 
geometrical  approximation . 

Without  the  addition  of  dynamic  step-size  adjustment,  the 
program  is  limited  to  applications  in  which  the  diode  pump  light 
illuminates  approximately  180  degrees  of  the  periphery  of  the  rod,  and 
the  total  pathlength  difference  from  the  diode  to  the  rod  surface 
between  a  ray  propagating  along  the  radius  of  the  system  and  a  ray 
propagating  to  the  extremes  of  the  rod  (near  where  rays  are  reflected  or 
refracted  such  that  they  miss  the  rod)  is  not  more  than  about  4000  or 
5000  wavelengths.  The  program  presumes  a  Gaussian  (in  one  dimension) 
far-field  pattern  from  the  diode. 

Computational  Example 

The  computation  is  illustrated  for  the  following  case: 


r i  =  0.3175  cm 

nl  = 

1.823 

r 2  =  0.5  cm 

n2  = 

1.333 

rj  =  0.65  cm 

n3  = 

1 . 5 

R  =  4  cm 

n  = 

1. 

The 

absorption  coefficient  in 

the  rod 

is  chosen  to 

be  a  =  3 

cm  1  ,  and 

the 

refectivity  at  the  wall 

is 

zero, 

because  the 

caustic  is  a  single 

pass 

phenomenon.  The  number 

of 

sources  is  chosen 

to  be  one,  because 

each 

source  produces  its  own  pair 

of  caustics .  The 

source  at 

R  =  4  (on 

0 

=  0)  produces  the  caustics 

running 

from  roughly 

(0.  235, ± 

140.85°)  to 

the 

rod  surface  at  (0.317 5, 

± 

164.03°).  About 

half-way 

along  the 

caustic  is  the  point  (0.27, 

± 

155.36 

°)  for  which 

we  made 

the  rough 

estimate  of  the  intensity  using  the  Airy  functions . 
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The  radii  and  refractive  indices  that  we  have  chosen  are  not 

from  a  point  design.  They  are  chosen  for  purposes  of  illustration 

only.  The  radiation  pattern  from  the  source  is  characterized  by  the 

2 

width  b Q  equal  to  6.436  degrees  in  the  formula  exp  -  (<5/$^)  ,  since 

this  particular  value  optimizes  the  flatness  of  the  pump  intensity  dis¬ 
tribution  when  multiple  pump  sources  are  used  with  a  wall  reflectivity 
of  0. 97. 


The  amplitude  and  phase  of  the  pump  light  just  inside  the 
surface  of  the  rod  on  the  side  facing  the  pump  source  are  readily 
tabluated  as  a  function  of  Q  using  the  two-dimensional  G.O.  theory,  the 
phase  being  simply  2ir/X  times  the  path  length  S2;  and  it  is  adjusted  to 
be  zero  at  the  entrance  point  on  the  central  ray  i.e.  at  Qq  in  Figure 
19.  The  path  length  3^(6)  is  as  shown  in  Fig.  20  and  the  amplitude  A(Q) 
is  shown  in  Fig.  21.  The  angle  9  for  the  ray  that  passes  through  the 
caustic  at  radius  0.27  cm  is  65.82  deg.  (1.15  radians) . 

The  range  of  0  in  Eq.  108a  is  limited  at  ±  93.76  degrees  where 
the  angle  of  incidence  8'  at  Q  reaches  90  degrees.  In  the  same 
equation,  the  obliquity  factor  has  to  take  into  account  the  fact  that 
the  surface  on  which  the  diffracting  field  is  defined  is  not  a  wave 
front;  and  we  do  this  by  using  the  standard  form  }  (1  +  cos  e)  for 
0(6)  and  multiplying  by  cos  8^  to  convert  the  surface  area 

element  to  the  wavefront  surface  area,  8^  being  the  angle  of  refraction . 

Figure  22a  shows  the  intensity  distribution  along  part  of  the 
circular  arc  a  =  0.27  computed  using  Eq.  108,  and  that  using  the  G.O. 
theory.  The  G.O.  curve  goes  to  infinity  at  about  155.36  deg.,  and 
elsewhere  approximates  the  mean  value  about  which  the  wave  optics 
intensity  oscillates.  The  peak  intensity  is  1.7634,  which  exceeds  the 
G.O.  value  for  the  same  point  by  0.7776,  and  which  is  a  factor  10.74 
greater  than  the  value  on  the  axis.  The  increment  0.7776  thus  deduced 
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Figure  21.  Plot  of  the  amplitude  of  the  source  wave  used  in  the 

diffraction  integral,  computed  just  inside  the  surface  of 
the  rod  (radius  0.3175  cm),  as  a  function  of  e  . 
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for  single  pass  (zero  wall  reflectivity)  can  now  be  added  as  a  cor¬ 
rection  to  the  G.O.  intensity  due  to  one  of  the  sources  for  multiple 
source  illumination  with  any  wall  reflectivity .  For  example,  with  five 
equally  spaced  sources  and  wall  reflectivity  equal  to  0.97,  the  G.O. 
illumination  at  the  point  (0.27,  155.445°)  is  1.640  and  at  (0,0)  is 

0.9388.  Adding  the  increment  to  the  first  number  and  dividing  by  the 
second  gives  a  ratio  of  2.58,  compared  to  the  10.63  deduced  for  a  single 
source,  i.e.  the  intensity  at  this  peak  is  2.58  times  the  intensity  on 
the  axis  of  the  rod.  (Of  course  this  result  depends  on  the  chosen  pump 
absorption  coefficient,  in  this  case  3  cm~^ . ) 

The  caustic  intensity  tends  to  increase  towards  the  edge  of 
the  rod,  this  tendency  being  offset  by  the  increased  attenuation  in  the 
medium.  Fig  22b  illustrates  the  computation,  again  for  a  single  pump 
source.  For  five  sources,  and  wall  relectivity  of  0.97,  the  enhancement 
over  the  intensity  on  the  axis  at  the  caustic  peak  near  the  edge  is  a 
factor  of  2.58. 

With  five  equally  spaced  sources  the  caustics  intersect  in 
planes  72  degrees  apart.  There  is  a  doubling  of  the  excess  intensity  at 
these  intersections .  Fig.  23  shows  a  radial  cut  through  one  caustic. 
The  excess  over  the  G.O.  intensity  is  doubled  and  added  to  the  G.O. 
result,  giving  in  this  case  an  enhancement  factor  of  3.80  over  the 
intensity  on  axis  with  wall  ref leu  ivity  0.97. 

The  Realistic  Case 


Our  model  computations  of  caustic  intensity  reveal  the 
expected  fringes  produced  by  interference  between  the  two  contributing 
optical  paths  The  fringe  spacing  becomes  smaller  as  one  goes  away  from 
the  G.O.  caustic,  as  the  angle  between  the  contributing  rays 
increases .  Because  the  fringe  spacing  and  the  G.O.  caustic  location  are 
both  sensitive  to  the  exact  geometrical  conditions  the  fringes  in  a  real 
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Intensity  vs.  Radius  at  0  -  144 
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Figure  23.  Radial  cut  of  intensity  profile  for  a  single  source.  The 
angle  144  degrees  would  be  at  a  crossing  of  caustics  from 
neighboring  sources  if  five  were  used. 


79 


system  will  tend  to  become  washed  out.  A  principal  contributing  factor 
is  the  longitudinal  propagation,  i.e.  the  fact  that  the  light  through  a 
field  point  originates  from  sources  distributed  along  z  emitting  over  a 
finite  angular  (6)  range. 

The  sensitivity  of  caustic  location  to  die  angle  (& )  is 
readily  demonstrated  as  follows.  For  the  standard  computational  example 
that  we  have  chosen  above  the  angular  location  of  the  G.O.  caustic  at  a 
=  0.27  cm  is  155 . 362288  deg.  This  corresponds  to  transverse  propaga¬ 
tion,  6=0.  When  6  has  some  other  value,  e.g.  1  deg.,  the  refractive 
indices  ( 1.823,  1.333,  1.5,  1)  are  replaced  according  to  Bq.  43  by  a  new 
set  (1.823194,  1.333089,  1.500127,  1);  and  with  these  new  parameters  the 
angular  position  of  the  G.O.  caustic  changes  to  155.369126  deg.  For  1 

deg.  steps  in  6  up  to  10  deg.  the  angular  position  of  the  caustic  a  at 

'  c 

radius  0.27  cm  is  found  to  satisfy  the  following  quadratic  expression  in 
6  : 

a  =  a  (o)  +  0.00685362 
c  c 

where  6  is  in  degrees.  For  example  5  =  3.82  deg.  will  move  the  caustic 
at  a  =  0.27  through  0.10  deg.  which  is  the  displacement  from  the  princi¬ 
pal  maximum  to  the  first  minimum.  Thus  6  ranges  of  a  few  degrees  are 
sufficient  to  wash  the  fringes  out,  producing  an  intensity  that  closely 
matches  the  value  predicted  by  G.O. 

The  important  problem  is  the  istimation  of  the  peak  inten¬ 
sity.  The  Airy  function  calculation  that  produced  Eq.  114  overestimated 
the  peak;  so  that  method  is  unsuitable.  The  theory  assumed  a  uniform 
amplitude  distribution  at  the  surface  of  the  rod,  instead  of  the  actual 
distribution  shown  in  Fig.  21.  Replacing  the  distribution  by  a  linear 
one  corresponding  to  the  tangent  to  the  amplitude  curve  will  also  give  a 
closed  form  mathematical  expression  for  the  intensity,  but  computational 
effort  is  still  needed  to  get  the  a  and  b  values  in  the  phase 
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expression.  ( See  Eq.  112.)  It  is  easier  to  run  PTE's  program  for  the 
integral  at  a  few  points  on  a  line  that  intersects  the  caustic,  and  get 
the  intensity  on  the  first  peak. 

Figs.  22  and  23  show  that  the  peak  intensity  is  of  the  order 
of  1.5  to  2  times  the  G.O.  value  at  the  same  point.  Thus  the  problem 
reduces  to  determining  the  distance  from  the  G.O.  caustic  (which  occurs 
at  the  inflection  point  on  the  Airy  Function  type  curve )  to  the  first 
peak.  Again  this  would  seem  to  require  determining  the  position  depen¬ 
dence  of  the  phase  coefficients  a  and  b.  We  know  of  no  short  cut,  like 
determining  a  Fresnel  number  for  example,  that  will  tell  use  the  posit¬ 
ion  at  which  to  compute  the  G.O.  intensity.  (The  Fresnel  number  comes 
out  of  a  quadratic  phase  dependence  rather  than  a  cubic  one.) 
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